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CONTROL SYSTEM AND METHOD, AND CONTROL UNIT 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention relates to a control system 
and method and an engine control unit that control a 
plant, using a model defining the relationship between 
a simulation value simulating the behavior of an 
internal variable of the plant and a detection value 
reflecting the behavior of the internal variable. 

Description of the Related Art 

Recently, due to social requirements, it is 
demanded of internal combustion engines that the 
engines have excellent exhaust emission characteristics, 
that is, an excellent emission reduction rate of the 
catalyst. On the other hand, internal combustion 
engines having a plurality of cylinders can suffer 
variation in air-fuel ratio between the cylinders to 
which the air- fuel mixture is supplied, due to the 
malfunction of an EGR system, an evaporative fuel 
processing system, or injectors. In such a case, there 
is a fear that the emission reduction rate of the 
catalyst is lowered. As a control system for a plant, 
which overcomes such a problem, there has been 
conventionally proposed an air-fuel ratio control 
system for an internal combustion engine, which 
corrects variation in air-fuel ratio between cylinders, 
using an observer based on the optimal control theory 
applied thereto (see e.g. Publication of Japanese 
Patent No. 3296472, pages 19-23, FIGS. 35 and 36). 
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This air-fuel ratio control system is comprised of a 
LAF sensor disposed in the collecting section of an 
exhaust pipe, for detecting the air-fuel ratio of 
exhaust gases, a control unit to which a detection 
signal (indicative of the detected air-fuel ratio) from 
the LAF sensor is input, and injectors disposed in the 
intake manifold of the exhaust pipe for the respective 
cylinders and connected to the control unit. 

In this control unit, variation in air-fuel ratio 
of exhaust gases emitted from a plurality of cylinders, 
i.e. variation in air-fuel ratio of the mixture between 
the cylinders is corrected by calculating a cylinder- 
try- cylinder fuel injection amount as the amount of fuel 
to be injected from each injector into the associated 
cylinder, based on the detected air-fuel ratio output 
from the LAF sensor, using the observer and by PID 
control, as described below. 

That is, the control unit calculates the basic 
fuel injection amount depending on the operating 
conditions of the engine, and multiplies the basic fuel 
injection amount by various correction coefficients to 
calculate the output fuel injection amount. Then, as 
described in detail hereinbelow, the observer 
calculates a cylinder-by-cylinder estimated air-fuel 
ratio, and a cylinder-by- cylinder estimated feedback 
correction coefficient is determined by PID control 
based on the estimated cylinder-by-cylinder air-fuel 
ratio. The cylinder-by-cylinder fuel injection amount 
is calculated by multiplying an output fuel injection 
amount by the cylinder-by- cylinder feedback correction 
coefficient . 

The cylinder-by-cylinder estimated air-fuel ratio 
is calculated by the observer based on the optimal 
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control theory. More specifically, by using a model of 
a discrete- time system representative of the 
relationship between a cylinder-by-cylinder fuel-air 
ratio and a fuel-air ratio detected at the collecting 
section (where the LAF sensor is disposed) , the 
cylinder-by-cylinder estimated air-fuel ratio is 
calculated. Further, in the PID control, a value 
obtained by dividing the fuel-air ratio detected at the 
collecting section, i.e. the detected air-fuel ratio, 
by the average value of the respective preceding values 
of the feedback correction coefficients is set to a 
target value, and the cylinder-by-cylinder feedback 
correction coefficient is calculated such that the 
difference between the target value and the cylinder- 
by-cylinder estimated air-fuel ratio calculated by the 
observer converges to a value of 0. 

Further, another air-fuel ratio control system is 
known which calculates the fuel injection amount on a 
cylinder-by-cylinder basis, based on an estimated 
intake air amount calculated by estimating the amount 
of intake air to be supplied to each of a plurality of 
cylinders, on a cylinder-by- cylinder basis, and an 
estimated air-fuel ratio calculated on a cylinder-by- 
cylinder basis by an observer similar to that described 
above (see e.g. Japanese Laid-Open Patent Publication 
(Kokai) No. 6-74076, pages, 3-12, FIGS. 1 and 31). 

More specifically, this air-fuel ratio control 
system calculates a target intake fuel amount by 
searching a map according to the engine speed and the 
intake pipe pressure. Further, by applying a fluid 
dynamics model to the intake system of the engine, the 
estimated intake air amount is calculated on a 
cylinder-by-cylinder basis, and the estimated air-fuel 
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ratio is calculated on a cylinder-by-cylinder basis, by 
the observer described above. Further, by dividing the 
estimated intake air amount by the estimated air- fuel 
ratio, an estimated intake fuel amount is calculated on 
a cylinder-by-cylinder basis, and a final fuel 
injection amount is calculated by an adaptive 
controller such that the estimated intake fuel amount 
becomes equal to the target intake fuel mount . 

Recently, aside from the above-mentioned demand 
of ensuring an excellent emission reduction rate of the 
catalyst, internal combustion engines are demanded of 
higher power output and higher torque. To meet the 
demand, there is conventionally employed the technique 
of reducing the exhaust resistance and exhaust 
interference by configuring the layout of the exhaust 
system into a complicated shape (in which exhaust 
passages from the cylinders are progressively combined 
in the exhaust manifold such that four passages, for 
example, are combined into two passages, and the two 
passages are then combined into one passage). However, 
when the conventional air-fuel ratio control system is 
applied to internal combustion engines having such a 
complicated exhaust system layout, the observer can no 
longer establish itself based on the conventional 
optimal control theory, and therefore, the variation in 
air-fuel ratio between the cylinders cannot be properly 
corrected, which can lead to a lowered emission 
reduction rate of the catalyst. This is because 
according to the conventional optimal control theory, 
modeling errors and changes in the dynamic 
characteristics of a model are not considered in the 
simulation model and the optimal control theory itself, 
which makes the observer small in margin of stability 
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and low in robustness. Therefore, the air-fuel ratio 
control system does not have a sufficient stability 
against changes in the contributions of exhaust gases 
from the individual cylinders to the detected air-fuel 
ratio of the LAF sensor caused by attachment of fuel, 
etc., changes in the response of the LAF sensor, and 
the aging of the same. 

Further, in the second-described air-fuel ratio 
control system, which uses the observer similar to that 
used in the first-described air-fuel ratio control 
system, there can be a case in which the observer 
cannot establish itself for the reason described above. 
In such a case, the fuel injection amount cannot be 
properly calculated on a cylinder-by- cylinder basis, 
which can degrade the emission reduction rate of the 
catalyst. Further, in a multi-cylinder internal 
combustion engine, in general, variation also occurs in 
intake air amount between the cylinders. However, the 
second-described air-fuel ratio control system does not 
consider the correction of the variation in intake air 
amount, and only estimates the intake air amount on a 
cylinder-by-cylinder basis, by applying the fluid 
dynamics model thereto. Therefore, the variation in 
intake air amount between the cylinders cannot be 
properly corrected, which brings about variation in the 
air- fuel ratio between the cylinders, causing further 
degradation of emission reduction rate of the catalyst. 



SUMMARY OF THE INVENTION 



It is a first object of the present invention to 
provide a control system and method and a control unit 
that are capable of realizing highly robust control 
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having a large margin of stability. 

It is a second object of the present invention to 
provide a control system and method and a control unit 
that are applicable to control of an air-fuel ratio of 
a mixture supplied to an internal combustion engine 
having a plurality of cylinders , and capable of 
appropriately and promptly correcting variation in air- 
fuel ratio or intake air amount between the cylinders 
and thereby realizing an accurate air-fuel ratio 
control even when the engine has a complicated exhaust 
system layout. 

To attain the first object, in a first aspect of 
the present invention, there is provided a control 
system for controlling a plant, comprising: 

detection means for detecting a detection value 
reflecting a behavior of a first internal variable of 
the plant; 

simulation value-generating means for generating 
a simulation value simulating the behavior of the first 
internal variable ; 

estimation means for estimating an estimation 
value of the detection value based on a model defining 
a relationship between the estimation value and the 
s imulat ion value ; 

identification means for identifying a model 
parameter of the model according to the detected 
detection value and the generated simulation value, 
such that the estimated estimation value becomes equal 
to the detected detection value; and 

first control means for determining a first input 
to be inputted to the plant, according to the 
identified model parameter. 

With the arrangement of the control system 
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according to the first aspect of the invention, the 
detection value reflecting the behavior of the first 
internal variable of the plant is detected, and the 
estimation value of the detection value is estimated 
based on a model defining the relationship between the 
estimation value and the simulation value simulating 
the behavior of the first internal variable. The model 
parameter of the model is identified according to the 
detection value and the simulation value, such that the 
estimated estimation value becomes equal to the 
detected detection value, and the first input to be 
inputted to the plant is determined according to the 
identified model parameter. Thus, the model parameter 
is identified such that the estimated estimation value 
becomes equal to the detected detection value, which 
enables the model parameter to be identified as a value 
in which the actual behavior of the first internal 
variable is properly reflected, particularly, enables 
the model parameter to be identified as a value in 
which the actual behavior of the first internal 
variable is reflected in real time, when an onboard 
identifier is used as the identification means. 
Further, the first input is determined according to the 
thus identified model parameter, so that even when the 
first internal variable is drastically changed, the 
first input can be determined as a value in which the 
behavior of the first internal variable is promptly and 
properly reflected, and by using the first input thus 
determined, it is possible to promptly and properly 
control the first internal variable to a predetermined 
state or a predetermined value. As a result, e.g. when 
the plant is to be controlled such that the first input 
causes the detection value detected by the detection 
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means to converge to a predetermined target value, even 
if the S/N ratio or sensitivity of the detection means 
is low, it is possible to set the detection value 
susceptible to the behavior of the first internal 
variable to the predetermined target value promptly 
with stability by causing the behavior of the first 
internal variable to be reflected in the first input. 
That is, it is possible to realize a control having a 
higher robustness and a larger margin of stability than 
the prior art . 

Preferably, the control system further comprises 
second control means for determining a second input to 
be inputted to the plant such that the detection value 
is caused to converge to a predetermined target value, 
the first internal variable comprising a plurality of 
first internal variables, the simulation value 
comprising a plurality of simulation values simulating 
respective behaviors of the plurality of first internal 
variables, the model parameter comprising a plurality 
of model parameters, and the identification means 
identifies the plurality of model parameters according 
to the detection value and the plurality of simulation 
values such that the estimated estimation value becomes 
equal to the detected detection value, the first 
control means determining the first input such that the 
identified model parameters converge to an average 
value thereof . 

With the arrangement of the preferred embodiment , 
the second control means determines the second input to 
be inputted to the plant such that the detection value 
is caused to converge to the predetermined target value, 
and the identification means identifies the plurality 
of model parameters according to the detection value 
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and the plurality of simulation values such that the 
estimated estimation value becomes equal to the 
detected detection value. The first control means 
determines the first input such that the identified 
model parameters converge to the average value thereof. 
Thus, the first input is determined such that the 
identified values of the plurality of model parameters 
converge to the average value thereof, which makes it 
possible to prevent a control process for causing the 
detection value detected by the detection means to 
converge to the predetermined target value and a 
control process for controlling the first internal 
variable from interfering with each other, and at the 
same time correct variation in behavior between the 
plurality of first internal variables. 

More preferably, the first control means 
comprises learned correction value-calculating means 
for calculating a learned correction value of the first 
input, using a sequential statistical algorithm, 
correction means for correcting the first input using 
the calculated learned correction value, and input 
means for inputting the corrected first input to the 
plant . 

The least -squares method is generally employed as 
the identifying computational algorithm. However, in 
the identifying computation by the least -squares method, 
after collecting a plurality of numbers of various data 
for computation, the computation is executed 
collectively based on the collected data. Therefore, 
at the start of the control, the identification of the 
model parameter is not executed until completion of 
collection of the data, which makes it impossible to 
calculate the first input based on the identified value 
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of the model parameter, which can degrade the 
controllability. In contrast, with the arrangement of 
the present preferred embodiment of the control system, 
the learned correction value of the first input is 
calculated with the sequential statistical algorithm, 
which enables the first input to be corrected even at 
the start of the control by the learned correction 
value calculated every control cycle. Therefore, e.g. 
by setting an initial value of the first input in 
advance, even before the model parameter is newly 
identified at the start of the control, the first input 
can be always corrected by the learned correction value 
calculated every control cycle, whereby the 
controllability at the start of the control can be 
enhanced. 

More preferably, the learned correction value - 
calculating means calculates the learned correction 
value of the first input using a regression equation in 
which the learned correction value is used as a 
dependent variable and a second internal variable 
having influence on the first internal variable is used 
as an independent variable, and calculates a regression 
coefficient and a constant term of the regression 
equation with the sequential statistical algorithm. 

With the arrangement of the preferred embodiment, 
the learned correction value of the first input is 
calculated using the regression equation in which the 
learned correction value is used as the dependent 
variable and a second internal variable having 
influence on the first internal variable is used as the 
independent variable, and the regression coefficient 
and the constant term of the regression equation are 
calculated with the sequential statistical algorithm. 
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Therefore, even when the rate of change in the second 
internal variable is very high, making the rate of 
change in the first internal variable also so high that 
it is difficult to estimate the first internal variable, 
it is possible to calculate the learned correction 
value as a value in which the actual state of the first 
internal variable is properly reflected, thereby 
further enhancing the controllability of the first 
internal variable by the first input. 

Preferably, the first control means determines an 
input component contained in the first input based on a 
difference between the model parameter and a 
predetermined target value. 

With the arrangement of this preferred embodiment, 
it is possible to determine the input component 
contained in the first input based on the difference 
between the model parameter and the predetermined 
target value. Therefore, the plant can be controlled 
such that model parameter converges to a predetermined 
target value, thereby causing the first internal 
variable of the plant to converge to a predetermined 
value without causing a steady- state deviation. 

More preferably, the first control means 
determines other input components than the input 
component contained in the first input, based on the 
model parameter. 

With the arrangement of this preferred embodiment, 
the first input contains not only the input component 
determined based on the difference between the model 
parameter and the predetermined target value, but also 
the other input components determined based on the 
model parameter. Therefore, e.g. when the plant is 
controlled such that the model parameter converges to 
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the predetermined target value, the first internal 
variable of the plant can be controlled that it 
converges to the predetermined value without causing 
overshooting or an oscillatory behavior. As a result, 
the detection value can be controlled to the stable 
state while preventing the same from becoming 
oscillatory or being overshot. 

Preferably, the first control means determines 
the first input according to the model parameter with a 
response-specified control algorithm. 

With the arrangement of this preferred embodiment, 
the first input is determined according to the model 
parameter with the response- specif ied control algorithm, 
and therefore, it is possible to control the plant, for 
example, such that model parameter converges to the 
predetermined target value, whereby the first internal 
variable of the plant can be controlled such that it 
converges to the predetermined value without causing 
overshooting or an oscillatory behavior. As a result, 
when the plant is controlled by the first input, the 
detection value can be controlled to a stable state 
while preventing the same from becoming oscillatory or 
overshot . 

Preferably, the identification means identifies 
the model parameter by a fixed gain method. 

With the arrangement of the preferred embodiment , 
the model parameter is identified by the fixed gain 
method, and therefore, it is possible to reduce 
computational load on the identification means. This 
makes it possible to shorten the computing time of the 
first input, whereby it is possible to calculate the 
first input promptly and properly as a value in which 
the behavior of the first internal variable is properly 
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reflected, even when the rate of change in the first 
internal variable is high. Further, when a method of 
identifying the model parameter by adding a 
predetermined correction component to the reference 
value thereof is employed as the fixed gain method, the 
identified value of the model parameter can be 
constrained to values close to the reference value, 
which makes it possible to prevent an increase in the 
rate of change in the first internal variable from 
causing the state of the first internal variable to be 
unsuitably reflected in the identified value of the 
model parameter, thereby making it possible to enhance 
the stability of the control. 

Further preferably, the identification means 
identifies the model parameter by calculating a model 
parameter reference value according to the second 
internal variable, and adding a predetermined 
correction component to the calculated model parameter 
reference value. 

With the arrangement of this preferred embodiment, 
the model parameter is identified by adding the 
predetermined correction component to the model 
parameter reference value calculated according to the 
second internal variable. This makes it possible to 
constrain the identified value of the model parameter 
to values close to the model parameter reference value, 
whereby even when the rate of change in the first 
internal variable is high due to the influence of 
change in the second internal variable, it is possible 
to promptly and properly calculate the first input as a 
value in which the behavior of the first internal 
variable is properly reflected, thereby enhancing the 
stability of the control. 
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Preferably, the control system further comprises 
delay means for delaying one of the detection value and 
the simulation value by a predetermined delay time 
period, and the identification means identifies the 
model parameter according to the delayed one of the 
detection value and the simulation value, and the other 
of the detection value and the simulation value. 

With the arrangement of this preferred embodiment, 
the model parameter is identified according to the 
delayed one of the detection value and the simulation 
value, and the other of the detection value and the 
simulation value. Therefore, e.g. when the detection 
value or the simulation value surfers from the dead 
time, it is possible to identify the model parameter 
with accuracy while taking the dead time into account, 
thereby further enhancing the stability of the control. 

Preferably the control system further comprises 
filter means for generating a filtered value of the 
detection value by subjecting the detection value to 
predetermined filtering processing, and the 
identification means identifies the model parameter 
according to the filtered value of the detection value 
and the simulation value. 

In general, in this kind of control system, when 
the absolute value of the detection value changes over 
a wide range, the identifying process by the 
identification means can be incapable of following up 
the change in the detection value, which can cause 
delay in identification of the model parameter, causing 
degraded accuracy of the identification. In contrast, 
with the arrangement of this preferred embodiment, the 
identification means identifies the model parameter 
according to the filtered value of the detection value 
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obtained by subjecting the detection value to the 
predetermined filtering processing and the simulation 
value, and therefore, by properly setting the filtering 
characteristics of the filtering processing, it is 
possible, even when the absolute value of the detection 
value changes over a wide range, the filtered value of 
the detection value can be generated as a value which 
positively contains information necessary for 
identification of the model parameter, i.e. information 
indicative of the behavior of the internal variables, 
and suppressed in the range of change thereof. 
Therefore, by identifying model parameter using the 
filtered value and the simulation value, it is possible 
to suppress delay in the identification of the model 
parameter and enhance the accuracy of the 

identification, thereby further enhancing the stability 
and response of the control. 

To attain the second object, in a second aspect 
of the present invention, there is provided a control 
system for an internal combustion engine including a 
plurality of cylinders, a plurality of exhaust passages 
extending from the plurality of cylinders, respectively, 
and one exhaust passage into which the plurality of 
exhaust passages are combined, the control system 
controlling an amount of fuel to be supplied to the 
plurality of cylinders, on a cylinder-by-cylinder basis, 
thereby controlling an air-fuel ratio of exhaust gases 
emitted from the plurality of cylinders, 

the control system comprising: 

fuel amount -determining means for determining an 
amount of fuel to be supplied to each of the plurality 
of cylinders; 

air-fuel ratio parameter-detecting means for 
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detecting an air-fuel ratio parameter indicative of an 
air- fuel ratio of exhaust gases in the one exhaust 
passage; 

simulation value-generating means for generating 
a plurality of simulation values simulating respective 
behaviors of air-fuel ratios of exhaust gases emitted 
from the plurality of cylinders ; 

estimation means for estimating an estimation 
value of the air-fuel ratio parameter based on a model 
defining a relationship between the estimation value 
and the plurality of simulation values; 

identification means for identifying a plurality 
of model parameters of the model according to the 
detected air-fuel ratio parameter and the generated 
plurality of simulation values, such that the 
estimation value of the air-fuel ratio parameter 
becomes equal to the detected air- fuel ratio parameter; 

first correction value-calculating means for 
calculating a first correction value for correcting the 
amount of fuel to be supplied to the plurality of 
cylinders, according to the identified plurality of 
model parameters, on a cylinder-by- cylinder basis; and 

first fuel amount -correcting means for correcting 
the determined amount of fuel according to the 
calculated first correction value, on a cylinder-by- 
cylinder basis . 

With the arrangement of the control system 
according to the second aspect of the invention, the 
amount of fuel to be supplied to each of the plurality 
of cylinders is determined by the fuel amount - 
determining means, and the air-fuel ratio parameter 
indicative of the air-fuel ratio of exhaust gases in 
the one exhaust passage is detected by the air-fuel 



H02-2749US 



17 



ratio parameter-detecting means. The estimation value 
of the air-fuel ratio parameter is estimated based on 
the model defining the relationship between the 
estimation value and the plurality of simulation values 
simulating respective behaviors of air-fuel ratios of 
exhaust gases emitted from the plurality of cylinders, 
and the plurality of model parameters of the model are 
identified by the identification means such that the 
estimation value of the air-fuel ratio parameter 
becomes equal to the detected air- fuel ratio parameter. 
The first correction value for correcting the amount of 
fuel to be supplied to the plurality of cylinders is 
calculated according to the identified plurality of 
model parameters, on a cylinder-by-cylinder basis, by 
the first correction value-calculating means. The 
determined fuel amount is corrected according to the 
calculated first correction value, on a cylinder-by- 
cylinder basis, by the first fuel amount -correcting 
means. Thus, the plurality of model parameters are 
identified such that the estimation value of the air- 
fuel ratio parameter becomes equal to the detected air- 
fuel ratio parameter, which makes it possible to 
identify the plurality of model parameters as values in 
which the actual behaviors of exhaust gases emitted 
from the plurality of cylinders, i.e. variation in air- 
fuel ratio between the cylinders is reflected therein. 
Therefore, by correcting the amount of fuel to be 
supplied to each cylinder according to the first 
correction value calculated according to the identified 
values of the plurality of model parameters, on a 
cylinder-by-cylinder basis, it is possible to properly 
correct variation in air fuel ratio between the 
cylinders. Further, by using an onboard identifier as 
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the identification means, it is possible to calculate 
the first correction value based on the model 
parameters identified in real time. This makes it 
possible, differently from the conventional control 
system, even when the dynamic characteristics of the 
controlled object are changed due to changes in 
respective contributions of the cylinders to the 
detected air-fuel ratio parameter, which are caused by 
attachment of fuel in the cylinders , variation in the 
response of the air-fuel ratio parameter-detecting 
means, and aging of the same, to correct the amount of 
fuel such that variation in air-fuel ratio between the 
cylinders is corrected (absorbed) while causing changes 
in the dynamic characteristics of the controlled object 
to be reflected in the model. As a result, even when 
the control system is applied to an internal combustion 
engine having a complicated exhaust system layout, it 
is possible to properly and promptly correct variation 
in air-fuel ratio between the cylinders, and thereby 
control the air-fuel ratio with accuracy. That is, it 
is possible to realize a highly robust air-fuel ratio 
control having a large margin of stability, and thereby, 
when a catalyst is provided in the exhaust passage, 
maintain an excellent emission reduction rate of the 
catalyst . 

Preferably, the control system further comprises 
second correction value-calculating means for 
calculating a second correction value for correcting 
the amount of fuel to be supplied to each cylinder, 
such that the air-fuel ratio parameter is caused to 
converge to a predetermined target value, and second 
fuel amount -correcting means for correcting the amount 
of fuel to be supplied to each cylinder according to 
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the calculated second correction value, and the first 
correction value-calculating means calculates the first 
correction value, on a cylinder-by-cylinder basis, such 
that the identified plurality of model parameters 
converge to an average value thereof. 

With the arrangement of this preferred embodiment, 
the second correction value-calculating means 
calculates the second correction value for correcting 
the amount of fuel to be supplied to each cylinder, 
such that the air-fuel ratio parameter is caused to 
converge to the predetermined target value, and the 
second fuel amount -correcting means corrects the amount 
of fuel to be supplied to each cylinder according to 
the calculated second correction value. Further, the 
first correction value-calculating means calculates the 
first correction value, on a cylinder-by-cylinder basis, 
such that the identified plurality of model parameters 
converge to an average value thereof. Thus, the first 
correction value is calculated such that the identified 
plurality of model parameters converge to an average 
value thereof, and therefore it is possible to correct 
variation in air-fuel ratio between the cylinders, 
whereby it is possible to prevent the control process 
for causing the air-fuel ratio parameter to converge to 
a predetermined target value and the control process 
for correcting variation in air-fuel ratio between the 
cylinders from interfering with each other, thereby 
ensuring stability of the air-fuel ratio control. 

Preferably, the control system further comprises 
learned correction value-calculating means for 
calculating a learned correction value of the first 
correction value with a sequential statistical 
algorithm, on a cylinder-by-cylinder basis, and the 
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first fuel amount -correcting means corrects the amount 
of fuel further according to the calculated learned 
correction value, on a cylinder-by-cylinder basis. 

As described hereinbefore, although the least- 
squares method is generally employed as the identifying 
computational algorithm, in the identifying computation 
by this method, after collecting a plurality of numbers 
of various data for computation, the computation is 
executed collectively based on the collected data. 
Therefore, at the start of the air-fuel ratio control, 
the identification of the model parameter is not 
executed until completion of collection of the data. 
This makes it impossible to calculate the first 
correction value based on the identified value of the 
model parameter, which can degrade the controllability 
of the air-fuel ratio control. In contrast, with the 
arrangement of the present preferred embodiment of the 
control system, the learned correction value of the 
first correction value is calculated with the 
sequential statistical algorithm, which enables the 
first correction value to be corrected by the learned 
correction value calculated every control cycle even at 
the start of the air-fuel ratio control. Therefore, by 
setting the initial value of the first correction value 
in advance, or by using the learned correction value 
calculated in the preceding operation of the engine as 
the initial value of the learned correction value of 
the current operation, it is possible, even before the 
identification of the model parameter is started at the 
start of the air-fuel ratio control, to always correct 
the first correction value by the learned correction 
value calculated every control cycle, whereby the 
controllability at the start of the air-fuel ratio 
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control can be enhanced. This makes it possible, when 
a catalyst is provided in the exhaust passage, to 
enhance the emission reduction rate of the catalyst at 
the start of the air-fuel ratio control. 

More preferably, the control system further 
comprises operating condition parameter-detecting means 
for detecting an operating condition parameter 
indicative of an operating condition of the engine, and 
the learned correction value-calculating means 
calculates the learned correction value using a 
regression equation in which the learned correction 
value is used as a dependent variable and the detected 
operating condition parameter is used as an independent 
variable, and calculates a regression coefficient and a 
constant term of the regression equation with the 
sequential statistical algorithm. 

With the arrangement of the preferred embodiment, 
the learned correction value of the first correction 
value is calculated using a regression equation in 
which the learned correction value of the first 
correction value is used as a dependent variable and 
the detected operating condition parameter is used as 
an independent variable, and the regression coefficient 
and the constant term of the regression equation are 
calculated with the sequential statistical algorithm. 
Therefore, even when the engine is in a drastically 
changing operating condition, such as a transient 
operating condition, causing a sudden change of the 
air-fuel . ratio, which makes it difficult to estimate 
the air-fuel ratio, it is possible to calculate the 
learned correction value as a value in which the actual 
state of the air-fuel ratio of each cylinder is 
properly reflected, thereby further enhancing the 
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controllability of the air-fuel ratio control. 

Preferably, the first correction value- 
calculating means calculates a correction value 
component contained in the first correction value based 
on a difference between the identified model parameters 
and a predetermined target value. 

With the arrangement of this preferred embodiment, 
the first correction value-calculating means calculates 
the correction value component contained in the first 
correction value based on the difference between the 
identified model parameters and the predetermined 
target value. Therefore, the amount of fuel can be 
corrected such that model parameters converge to the 
predetermined target value, thereby providing control 
on the air-fuel ratio, on a cylinder-by-cylinder basis, 
such that the air-fuel ratio converges to a 
predetermined value without causing a steady- state 
deviation . 

More preferably, the first correction value- 
calculating means calculates other correction value 
components than the correction value component 
contained in the first correction value, based on the 
identified model parameters. 

With the arrangement of this preferred embodiment, 
the first correction value contains not only the 
correction value component determined based on the 
difference between the model parameters and the 
predetermined target value, but also the other 
correction value components determined based on the 
model parameters. Therefore, e.g. when the amount of 
fuel is controlled, on a cylinder-by-cylinder basis, 
such that the model parameters converge to the 
predetermined target value, the air-fuel ratio can be 
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controlled, on a cylinder-by-cylinder basis, such that 
it converges to the predetermined value without causing 
overshooting or an oscillatory behavior, with stability. 

Preferably, the first correction value- 
calculating means calculates the first correction value 
according to the model parameters with a response- 
specified control algorithm. 

With the arrangement of this preferred embodiment, 
the first correction value is determined according to 
the model parameters with the response-specified 
control algorithm, and therefore, it is possible to 
correct the amount of fuel, for example, such that 
model parameters converge to the predetermined target 
value, whereby the air-fuel ratio can be corrected, on 
a cylinder-by-cylinder basis, such that it converges to 
the predetermined value without causing overshooting or 
an oscillatory behavior, with stability. 

Preferably, the identification means identifies 
the model parameters by a fixed gain method. 

With the arrangement of the preferred embodiment, 
the model parameters are identified by the fixed gain 
method, and therefore, it is possible to reduce 
computational load on the identification means. This 
makes it possible to shorten the computing time of the 
first correction value, whereby it is possible to 
calculate the first correction value promptly and 
properly, on a cylinder-by- cylinder basis, as a value 
in which the behavior of the air-fuel ratio is properly 
reflected, even when the rate of change in the air- fuel 
ratio of each cylinder is high due to a transient 
operating condition of the engine. Further, when a 
method of identifying the model parameters by adding 
respective predetermined correction components to 
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reference values thereof is employed as the fixed gain 
method, the identified values of the model parameters 
can be constrained to values close to the reference 
values, which makes it possible to prevent an increase 
in the rate of change in the air- fuel ratio from 
causing the actual state of the air-fuel ratio to be 
unsuitably reflected in the identified values of the 
model parameters, thereby making it possible to further 
enhance the stability of the air-fuel ratio control. 

Preferably, the identification means identifies 
the model parameters by calculating respective model 
parameter reference values according to the operating 
condition parameter, and adding predetermined 
correction components to the calculated model parameter 
reference values , respectively. 

With the arrangement of this preferred embodiment, 
the model parameters are identified by adding the 
respective predetermined correction components to the 
model parameter reference values calculated according 
to the operating condition parameter. This makes it 
possible to constrain the identified values of the 
model parameters to respective values close to the 
model parameter reference values, whereby even when the 
rate of change in the air-fuel ratio is high due to the 
influence of change in the operating condition of the 
engine, it is possible to promptly and properly 
calculate the first correction value, on a cylinder-by- 
cylinder basis, as a value in which the behavior of the 
air-fuel ratio is properly reflected, thereby further 
enhancing the stability of the control. 

Preferably, the control system further comprises 
delay means for delaying the air-fuel ratio parameter 
by a predetermined delay time period, and the 
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identification means identifies the model parameters 
according to the delayed air-fuel ratio parameter and 
the plurality of simulation values. 

In general, in the internal combustion engine, 
there is a predetermined dead time from a time that the 
mixture supplied to each cylinder has burned to a time 
that the resulting exhaust gases reach the collecting 
section of the exhaust passage or a location downstream 
of the collecting section. However, with the 
arrangement of this preferred embodiment, the model 
parameter is identified according to the delayed air- 
fuel ratio parameter, which is delayed by the 
predetermined dead time, and the plurality of model 
parameters. Therefore, it is possible to identify the 
model parameter with accuracy while taking the dead 
time into account, thereby further enhancing the 
stability of the control. 

To attain the second object, in a third aspect of 
the present invention, there is provided a control 
system for an internal combustion engine including one 
intake passage, a plurality of intake passages 
branching from the one intake passage, and a plurality 
of cylinders connected to the plurality of intake 
passages extend, respectively, the control system 
controlling an amount of fuel to be supplied to the 
plurality of cylinders, on a cylinder-by-cylinder basis, 
thereby controlling an air-fuel ratio of exhaust gases 
emitted from the plurality of cylinders, 

the control system comprising: 

fuel amount -determining means for determining an 
amount of fuel to be supplied to each of the plurality 
of cylinders; 

intake air amount parameter-detecting means 
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disposed in the one intake passage, for detecting an 
intake air amount parameter indicative of an amount of 
intake air; 

simulation value-generating means for generating 
a plurality of simulation values simulating respective 
behaviors of amounts of intake air to be drawn into the 
plurality of cylinders ; 

estimation means for estimating an estimation 
value of the intake air amount parameter based on a 
model defining a relationship between the estimation 
value and the plurality of simulation values; 

identification means for identifying a plurality 
of model parameters of the model according to the 
detected intake air amount parameter and the generated 
plurality of simulation values, such that the 
estimation value of the intake air amount parameter 
becomes equal to the detected intake air amount 
parameter; 

third correction value-calculating means for 
calculating a third correction value for correcting the 
amount of fuel to be supplied to the plurality of 
cylinders, according to the identified plurality of 
model parameters, on a cylinder-by-cylinder basis; and 

third fuel amount -correcting means for correcting 
the determined amount of fuel according to the 
calculated third correction value, on a cylinder-by- 
cylinder basis. 

With the arrangement of the control system 
according to the third aspect of the invention, the 
fuel amount -determining means determines the amount of 
fuel to be supplied to each cylinder, and the intake 
air amount parameter-detecting means disposed in the 
one intake passage detects the intake air amount 
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parameter indicative of the amount of intake air. The 
estimation means estimates the estimation value of the 
intake air amount parameter based on the model defining 
the relationship between the estimation value and the 
plurality of simulation values simulating respective 
behaviors of amounts of intake air to be drawn into the 
plurality of cylinders, and the identification means 
identifies the plurality of model parameters of the 
model such that the estimation value of the intake air 
amount parameter becomes equal to the detected intake 
air amount parameter. The third correction value- 
calculating means calculates the third correction value 
for correcting the amount of fuel to be supplied to the 
plurality of cylinders, according to the identified 
plurality of model parameters, on a cylinder-by- 
cylinder basis, and the third fuel amount -correcting 
means corrects the determined fuel amount according to 
the calculated third correction value, on a cylinder- 
by-cylinder basis. Thus, the plurality of model 
parameters are identified such that the estimation 
value of the intake air amount parameter becomes equal 
to the detected intake air amount parameter, which 
makes it possible to identify the plurality of model 
parameters as values in which the actual behaviors of 
amounts of intake air drawn into the cylinders are 
reflected therein, i.e. variation in intake air amount 
between the cylinders is reflected therein. Therefore, 
by correcting the amount of fuel to be supplied to each 
cylinder according to the third correction value 
calculated according to the identified values of the 
plurality of model parameters, on a cylinder-by- 
cylinder basis, it is possible to properly correct 
variation in intake air amount between the cylinders. 
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Further, by using an onboard identifier as the 
identification means, it is possible to calculate the 
third correction value based on the model parameters 
identified in real time. This makes it possible, 
differently from the conventional control system, even 
when the dynamic characteristics of the controlled 
object are changed due to variation in the response of 
the intake air amount parameter-detecting means, and 
aging of the same, to correct the fuel amount such that 
variation in intake air amount between the cylinders is 
corrected while causing changes in the dynamic 
characteristics of the controlled object to be 
reflected in the model. As a result, even when the 
control system is applied to an internal combustion 
engine having a complicated exhaust system layout, it 
is possible to properly and promptly correct variation 
in intake air amount between the cylinders, and thereby 
control the air-fuel ratio with accuracy. That is, it 
is possible to realize a highly robust air-fuel ratio 
control having a large margin of stability, and thereby, 
when a catalyst is provided in the exhaust passage, 
maintain an excellent emission reduction rate of the 
catalyst . 

Preferably, the engine includes a plurality of 
exhaust passages extending from the plurality of 
cylinders, respectively, and one exhaust passage into 
which the plurality of exhaust passages are combined, 
and the control system further comprises intake air 
amount parameter-detecting means for detecting an 
intake air amount parameter indicative of an air- fuel 
ratio of exhaust gases in the one exhaust passage, 
fourth correction value-calculating means for 
calculating a fourth correction value for correcting 
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the amount of fuel to be supplied to each cylinder, 
such that the detected air-fuel ratio parameter is 
caused to converge to a predetermined target value, and 
fourth fuel amount -correcting means for correcting the 
amount of fuel to be supplied to each cylinder 
according to the calculated fourth correction value, 
the third correction value-calculating means 
calculating the third correction value, on a cylinder- 
by-cylinder basis, such that the identified plurality 
of model parameters converge to an average value 
thereof . 

With the arrangement of this preferred embodiment, 
the fourth correction value-calculating means 
calculates the fourth correction value for correcting 
the amount of fuel to be supplied to each cylinder, 
such that the air-fuel ratio parameter is caused to 
converge to the predetermined target value, and the 
fourth fuel amount -correcting means corrects the amount 
of fuel to be supplied to each cylinder according to 
the calculated fourth correction value. Further, the 
third correction value-calculating means calculates the 
third correction value, on a cylinder-by-cylinder basis, 
such that the identified plurality of model parameters 
converge to an average value thereof . Thus , the third 
correction value is calculated such that the identified 
plurality of model parameters converge to an average 
value thereof, which makes it possible to correct 
variation in intake air amount between the cylinders, 
whereby it is possible to prevent the control process 
for causing the air-fuel ratio parameter to converge to 
the predetermined target value and the control process 
for correcting variation in intake air amount between 
the cylinders from interfering with each other, thereby 
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ensuring stability of the air-fuel ratio control. 

Preferably, the control system further comprises 
learned correction value-calculating means for 
calculating a learned correction value of the third 
correction value with a sequential statistical 
algorithm, on a cylinder-by-cylinder basis, and the 
third fuel amount -correcting means corrects the amount 
of fuel further according to the calculated learned 
correction value, on a cylinder-by-cylinder basis. 

As described hereinbefore, when the least -squares 
method is employed as the identifying computational 
algorithm, the identification of the model parameter is 
not executed until completion of collection of the data 
at the start of the control, which makes it impossible 
to calculate the third correction value based on the 
identified value of the model parameter, which can 
degrade the controllability of the air-fuel ratio 
control. In contrast, with the arrangement of the 
present preferred embodiment of the control system, the 
learned correction value of the third correction value 
is calculated with the sequential statistical algorithm, 
which enables the third correction value to be 
corrected by the learned correction value calculated 
every control cycle even at the start of the control. 
Therefore, by setting the initial value of the third 
correction value in advance, or by using the learned 
correction value calculated in the preceding operation 
of the engine as the initial value of the learned 
correction value of the current operation, it is 
possible, even before the identification of the model 
parameter is started at the start of the air-fuel ratio 
control, to always correct the third correction value 
by the learned correction value calculated every 
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control cycle, whereby the controllability at the start 
of the air-fuel ratio control can be enhanced. This 
makes it possible, when a catalyst is provided in the 
exhaust passage, to enhance the emission reduction rate 
of the catalyst at the start of the air-fuel ratio 
control . 

More preferably, the control system further 
comprises operating condition parameter-detecting means 
for detecting an operating condition parameter 
indicative of an operating condition of the engine, and 
the learned correction value-calculating means 
calculates the learned correction value using a 
regression equation in which the learned correction 
value is used as a dependent variable and the detected 
operating condition parameter is used as an independent 
variable, and calculates a regression coefficient and a 
constant term of the regression equation with the 
sequential statistical algorithm. 

With the arrangement of the preferred embodiment, 
the learned correction value of the third correction 
value is calculated using the regression equation in 
which the learned correction value is used as the 
dependent variable and the detected operating condition 
parameter is used as the independent variable, and the 
regression coefficient and the constant term of the 
regression equation are calculated with the sequential 
statistical algorithm. Therefore, even when the engine 
is in a drastically changing operating condition, such 
as a transient operating condition, causing a sudden 
change of the air- fuel ratio, which makes it difficult 
to estimate the first internal variable, it is possible 
to calculate the learned correction value as a value in 
which the actual state of the amount of intake air 
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supplied to each cylinder is properly reflected, 
thereby further enhancing the controllability of the 
air-fuel ratio control. 

Preferably, the third correction value- 
calculating means calculates a correction value 
component contained in the third correction value based 
on a difference between the identified model parameters 
and a predetermined target value. 

With the arrangement of this preferred embodiment , 
the third correction value-calculating means calculates 
the correction value component contained in the third 
correction value based on the difference between the 
identified model parameters and the predetermined 
target value. Therefore, the amount of fuel can be 
corrected such that the model parameters converge to a 
predetermined target value, thereby providing control 
on the intake air amount on a cylinder-by-cylinder 
basis such that the intake air amount converges to a 
predetermined value without causing a steady-state 
deviation . 

More preferably, the third correction value- 
calculating means calculates other correction value 
components than the correction value component 
contained in the third correction value, based on the 
identified model parameters. 

With the arrangement of this preferred embodiment, 
the third correction value contains not only the 
correction value component determined based on the 
difference between the model parameters and the 
predetermined target value, but also other correction 
value components determined based on the model 
parameters. Therefore, e.g. when the amount of fuel is 
controlled on a cylinder-by- cylinder basis such that 
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the model parameters converge to the predetermined 
target value , the amount of intake air can be 
controlled, on a cylinder-by- cylinder basis, such that 
it converges to the predetermined value without causing 
overshooting or an oscillatory behavior, with stability. 

Preferably, the third correction value - 
calculating means calculates the third correction value 
according to the model parameters with a response- 
specified control algorithm. 

With the arrangement of this preferred embodiment , 
the third correction value is determined according to 
the model parameters with the response-specified 
control algorithm, and therefore, it is possible to 
correct the amount of fuel, for example, such that 
model parameters converge to the predetermined target 
value, whereby the amount of intake air can be 
corrected, on a cylinder-by-cylinder basis, such that 
it converges to the predetermined value without causing 
overshooting or an oscillatory behavior, with stability. 

Preferably, the identification means identifies 
the model parameters by a fixed gain method. 

With the arrangement of the preferred embodiment, 
the model parameters are identified by the fixed gain 
method, and therefore, it is possible to reduce 
computational load on the identification means. This 
makes it possible to shorten the computing time of the 
third correction value, whereby it is possible to 
calculate the third correction value promptly and 
properly, on a cylinder-by-cylinder basis, as a value 
in which the behavior of the amount of intake air is 
properly reflected, even when the rate of change in the 
amount of intake air of each cylinder is high due to a 
transient operating condition of the engine. Further, 
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when a method of identifying the model parameters by 
adding respective predetermined correction components 
to reference values thereof is employed as the fixed 
gain method, the identified values of the model 
parameters can be constrained to values close to the 
reference values, which makes it possible to prevent an 
increase in the rate of change in the intake air amount 
from causing the actual state of the intake air amount 
to be unsuitably reflected in the identified values of 
the model parameters, thereby making it possible to 
further enhance the stability of the air-fuel ratio 
control . 

Further preferably, the identification means 
identifies the model parameters by calculating 
respective model parameter reference values according 
to the operating condition parameter, and adding 
predetermined correction components to the calculated 
model parameter reference values, respectively. 

With the arrangement of this preferred embodiment, 
the model parameters are identified by adding the 
respective predetermined correction components to the 
model parameter reference values calculated according 
to the operating condition parameter. This makes it 
possible to constrain the identified values of the 
model parameters to values close to the model parameter 
reference values, whereby even when the rate of change 
in the amount of intake air is high due to the 
influence of change in the operating condition of the 
engine, it is possible to promptly and properly 
calculate the third correction value, on a cylinder-by- 
cylinder basis, as a value in which the behavior of the 
amount of intake air is properly reflected, thereby 
further enhancing the stability of the control. 
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Preferably, the control system further comprises 
delay means for delaying the plurality of simulation 
values by a predetermined delay time period, and the 
identification means identifies the model parameters 
according to the intake air amount parameter and the 
delayed plurality of simulation values. 

In general, in the internal combustion engine, 
there is a predetermined dead time before air drawn 
into the intake passage reaches each cylinder via the 
branches of the intake passage. However, with the 
arrangement of this preferred embodiment, the model 
parameters are identified according to the intake air 
amount parameter and the plurality of delayed 
simulation values, which is delayed by the 
predetermined dead time. Therefore, it is possible to 
identify the model parameters with accuracy while 
taking the dead time into account, thereby further 
enhancing the stability of the control. 

Preferably, the control system further comprises 
filter means for generating a filtered value of the 
intake air amount parameter by subjecting the intake 
air amount parameter to predetermined filtering 
processing, and the identification means identifies the 
model parameters according to the generated filtered 
value of the intake air amount parameter and the 
plurality of simulation values. 

In general, in this kind of control system, when 
the engine is in an operating condition in which the 
absolute value of the intake air amount parameter 
changes over a wide range, such as a transient 
operating condition, the identifying process by the 
identification means can be incapable of following up 
the change, which can cause delay in identification of 
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the model parameters, causing degraded accuracy of the 
identification. In contrast, with the arrangement of 
this preferred embodiment, the identification means 
identifies the model parameters according to the 
filtered value of the intake air amount parameter 
obtained by subjecting the intake air amount parameter 
to predetermined filtering processing and the 
simulation values, and therefore, by properly setting 
the filtering characteristics of the filtering 
processing, it is possible, even when the absolute 
value of the intake air amount parameter changes over a 
wide range, the filtered value of the intake air amount 
parameter value can be generated as a value which 
positively contains information necessary for 
identification of the model parameter, i.e. information 
indicative of the behavior (variation and the like) of 
the intake air of each cylinder, and is suppressed in 
the range of change thereof. Therefore, by identifying 
using the filtered value and the simulation values, it 
is possible to suppress delay in the identification of 
the model parameters and enhance the accuracy of the 
identification, thereby further enhancing the stability 
and response of the air-fuel ratio control. 

To attain the first object, in a fourth aspect of 
the present invention, there is provided a control 
method for controlling a plant, comprising: 

a detection step of detecting a detection value 
reflecting a behavior of a first internal variable of 
the plant; 

a simulation value-generating step of generating 
a simulation value simulating the behavior of the first 
internal variable; 

an estimation step of estimating an estimation 



H02-2749US 

37 

value of the detection value based on a model defining 
a relationship between the estimation value and the 
simulation value; 

an identification step of identifying a model 
parameter of the model according to the detected 
detection value and the generated simulation value, 
such that the estimated estimation value becomes equal 
to the detected detection value; and 

a first control step of determining a first input 
to be inputted to the plant, according to the 
identified model parameter. 

With the arrangement of the control method 
according to the fourth aspect of the present invention, 
it is possible to obtain the same advantageous effects 
as provided by the first aspect of the present 
invention . 

Preferably, the control method further comprises 
a second control step of determining a second input to 
be inputted to the plant such that the detection value 
is caused to converge to a predetermined target value, 
the first internal variable comprising a plurality of 
first internal variables, the simulation value 
comprising a plurality of simulation values simulating 
respective behaviors of the plurality of first internal 
variables, the model parameter comprising a plurality 
of model parameters, and the identification step 
includes identifying the plurality of model parameters 
according to the detection value and the plurality of 
simulation values such that the estimated estimation 
value becomes equal to the detected detection value, 
the first control step including determining the first 
input such that the identified model parameters 
converge to an average value thereof . 
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Preferably, the first control step comprises a 
learned correction value-calculating step of 
calculating a learned correction value of the first 
input, using a sequential statistical algorithm, a 
correction step of correcting the first input using the 
calculated learned correction value, and an input step 
of inputting the corrected first input to the plant. 

More preferably, the learned correction value- 
calculating step includes calculating the learned 
correction value of the first input using a regression 
equation in which the learned correction value is used 
as a dependent variable and a second internal variable 
having influence on the first internal variable is used 
as an independent variable, and calculating a 
regression coefficient and a constant term of the 
regression equation with the sequential statistical 
algorithm. 

Preferably, the first control step includes 
determining an input component contained in the first 
input based on a difference between the model parameter 
and a predetermined target value. 

More preferably, the first control step includes 
determining other input components than the input 
component contained in the first input, based on the 
model parameter. 

Preferably, the first control step includes 
determining the first input according to the model 
parameter with a response-specified control algorithm. 

Preferably, the identification step includes 
identifying the model parameter by a fixed gain method. 

Further preferably, the identification step 
includes identifying the model parameter by calculating 
a model parameter reference value according to the 
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second internal variable, and adding a predetermined 
correction component to the calculated model parameter 
reference value. 

Preferably, the control method further comprises 
a delay step of delaying one of the detection value and 
the simulation value by a predetermined delay time 
period, and the identification step includes 
identifying the model parameter according to the 
delayed one of the detection value and the simulation 
value, and the other of the detection value and the 
simulation value. 

Preferably, the control method further comprises 
a filter step of generating a filtered value of the 
detection value by subjecting the detection value to 
predetermined filtering processing, and the 
identification step includes identifying the model 
parameter according to the filtered value of the 
detection value and the simulation value. 

With the arrangements of these preferred 
embodiments, it is possible to obtain the same 
advantageous effects as provided by the corresponding 
preferred embodiments of the first aspect of the 
present invention . 

To attain the second object, in a fifth aspect of 
the present invention, there is provided a control 
method for an internal combustion engine including a 
plurality of cylinders, a plurality of exhaust passages 
extending from the plurality of cylinders, respectively, 
and one exhaust passage into which the plurality of 
exhaust passages are combined, the control method 
controlling an amount of fuel to be supplied to the 
plurality of cylinders, on a cylinder-by-cylinder basis, 
thereby controlling an air-fuel ratio of exhaust gases 
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emitted from the plurality of cylinders, 
the control method comprising: 

a fuel amount -determining step of determining an 
amount of fuel to be supplied to, each of the plurality 
of cylinders; 

an air-fuel ratio parameter-detecting step of 
detecting an air-fuel ratio parameter indicative of an 
air-fuel ratio of exhaust gases in the one exhaust 
passage ; 

a simulation value-generating step of generating 
a plurality of simulation values simulating respective 
behaviors of air-fuel ratios of exhaust gases emitted 
from the plurality of cylinders; 

an estimation step of estimating an estimation 
value of the air-fuel ratio parameter based on a model 
defining a relationship between the estimation value 
and the plurality of simulation values; 

an identification step of identifying a plurality 
of model parameters of the model according to the 
detected air-fuel ratio parameter and the generated 
plurality of simulation values, such that the 
estimation value of the air- fuel ratio parameter 
becomes equal to the detected air-fuel ratio parameter; 

a first correction value-calculating step of 
calculating a first correction value for correcting the 
amount of fuel to be supplied to the plurality of 
cylinders , according to the identified plurality of 
model parameters, on a cylinder-by-cylinder basis; and 

a first fuel amount-correcting step of correcting 
the determined amount of fuel according to the 
calculated first correction value, on a cylinder-by- 
cylinder basis. 

With the arrangement of the control method 
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according to the fifth aspect of the present invention, 
it is possible to obtain the same advantageous effects 
as provided by the second aspect of the present 
invention . 

Preferably, the control method further comprises 
a second correction value-calculating step of 
calculating a second correction value for correcting 
the amount of fuel to be supplied to each cylinder, 
such that the air-fuel ratio parameter is caused to 
converge to a predetermined target value, and a second 
fuel amount -correcting step of correcting the amount of 
fuel to be supplied to each cylinder according to the 
calculated second correction value, and the first 
correction value-calculating step includes calculating 
the first correction value, on a cylinder-by-cylinder 
basis, such that the identified plurality of model 
parameters converge to an average value thereof . 

Preferably, the control method further comprises 
a learned correction value-calculating step of 
calculating a learned correction value of the first 
correction value with a sequential statistical 
algorithm,, on a cylinder-by-cylinder basis, and the 
first fuel amount -correcting step includes correcting 
the amount of fuel further according to the calculated 
learned correction value, on a cylinder-by-cylinder 
basis . 

More preferably, the control method further 
comprises an operating condition parameter-detecting 
step of detecting an operating condition parameter 
indicative of an operating condition of the engine, and 
the learned correction value-calculating step includes 
calculating the learned correction value using a 
regression equation in which the learned correction 
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value is used as a dependent variable and the detected 
operating condition parameter is used as an independent 
variable, and calculating a regression coefficient and 
a constant term of the regression equation with the 
sequential statistical algorithm. 

Preferably, the first correction value- 
calculating step includes calculating a correction 
value component contained in the first correction value 
based on a difference between the identified model 
parameters and a predetermined target value . 

More preferably, the first correction value- 
calculating step includes calculating other correction 
value components than the correction value component 
contained in the first correction value, based on the 
identified model parameters. 

Preferably, the first correction value- 
calculating step includes calculating the first 

a 

correction value according to the model parameters with 
a response- specif ied control algorithm. 

Preferably, the identification step includes 
identifying the model parameters by a fixed gain method. 

Further preferably, the identification step 
includes identifying the model parameters by 
calculating respective model parameter reference values 
according to the operating condition parameter, and 
adding predetermined correction components to the 
calculated model parameter reference values, 
respectively . 

Preferably, the control method further comprises 
a delay step of delaying the air-fuel ratio parameter 
by a predetermined delay time period, and the 
identification step includes identifying the model 
parameters according to the delayed air- fuel ratio 
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parameter and the plurality of simulation values. 

With the arrangements of these preferred 
embodiments, it is possible to obtain the same 
advantageous effects as provided by the corresponding 
preferred embodiments of the second aspect of the 
present invention . 

To attain the second object, in a sixth aspect of 
the present invention, there is provided a control 
method for an internal combustion engine including one 
intake passage, a plurality of intake passages 
branching from the one intake passage, and a plurality 
of cylinders connected to the plurality of intake 
passages extend, respectively, the control method 
controlling an amount of fuel to be supplied to the 
plurality of cylinders, on a cylinder-by-cylinder basis, 
thereby controlling an air-fuel ratio of exhaust gases 
emitted from the plurality of cylinders, 

the control method comprising: 

a fuel amount -determining step of determining an 
amount of fuel to be supplied to each of the plurality 
of cylinders; 

an intake air amount parameter-detecting step of 
detecting an intake air amount parameter indicative of 
an amount of intake air in the one intake passage; 

a simulation value-generating step of generating 
a plurality of simulation values simulating respective 
behaviors of amounts of intake air to be drawn into the 
plurality of cylinders; 

an estimation step of estimating an estimation 
value of the intake air amount parameter based on a 
model defining a relationship between the estimation 
value and the plurality of simulation values; 

an identification step of identifying a plurality 
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of model parameters of the model according to the 
detected intake air amount parameter and the generated 
plurality of simulation values, such that the 
estimation value of the intake air amount parameter 
becomes equal to the detected intake air amount 
parameter; 

a third correction value-calculating step of 
calculating a third correction value for correcting the 
amount of fuel to be supplied to the plurality of 
cylinders, according to the identified plurality of 
model parameters, on a cylinder-by-cylinder basis; and 

a third fuel amount -correcting step of correcting 
the determined amount of fuel according to the 
calculated third correction value, on a cylinder-by- 
cylinder basis. 

With the arrangement of the control method 
according to the sixth aspect of the present invention, 
it is possible to obtain the same advantageous effects 
as provided by the first aspect of the present 
invention . 

Preferably, the engine includes a plurality of 
exhaust passages extending from the plurality of 
cylinders, respectively, and one exhaust passage into 
which the plurality of exhaust passages are combined, 
and the control method further comprises an air-fuel 
ratio parameter-detecting step of detecting an air-fuel 
ratio parameter indicative of an air-fuel ratio of 
exhaust gases in the one exhaust passage, a fourth 
correction value-calculating step of calculating a 
fourth correction value for correcting the amount of 
fuel to be supplied to each cylinder, such that the 
detected air-fuel ratio parameter is caused to converge 
to a predetermined target value, and a fourth fuel 
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amount-correcting step of correcting the amount of fuel 
to be supplied to each cylinder according to the 
calculated fourth correction value, the third 
correction value-calculating step including calculating 
the third correction value, on a cylinder-by-cylinder 
basis, such that the identified plurality of model 
parameters converge to an average value thereof. 

Preferably, the control method further comprises 
a learned correction value-calculating step of 
calculating a learned correction value of the third 
correction value with a sequential statistical 
algorithm, on a cylinder-by- cylinder basis, and the 
third fuel amount -correcting step includes correcting 
the amount of fuel further according to the calculated 
learned correction value, on a cylinder-by-cylinder 
basis . 

More preferably, the control method further 
comprises an operating condition parameter-detecting 
step of detecting an operating condition parameter 
indicative of an operating condition of the engine, and 
the learned correction value -calculating step includes 
calculating the learned correction value using a 
regression equation in which the learned correction 
value is used as a dependent variable and the detected 
operating condition parameter is used as an independent 
variable, and calculating a regression coefficient and 
a constant term of the regression equation with the 
sequential statistical algorithm. 

Preferably, the third correction value - 
calculating step includes calculating a correction 
value component contained in the third correction value 
based on a difference between the identified model 
parameters and a predetermined target value. 
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More preferably, the third correction value - 
calculating step includes calculating other correction 
value components than the correction value component 
contained in the third correction value, based on the 
identified model parameters. 

Preferably, the third correction value- 
calculating step includes calculating the third 
correction value according to the model parameters with 
a response-specified control algorithm. 

Preferably, the identification step includes 
identifying the model parameters by a fixed gain method. 

Further preferably, the identification step 
includes identifying the model parameters by 
calculating respective model parameter reference values 
according to the operating condition parameter, and 
adding predetermined correction components to the 
calculated model parameter reference values, 
respectively . 

Preferably, the control method further comprises 
a delay step of delaying the plurality of simulation 
values by a predetermined delay time period, and the 
identification step includes identifying the model 
parameters according to the intake air amount parameter 
and the delayed plurality of simulation values . 

Preferably, the control method further comprises 
a filter step of generating a filtered value of the 
intake air amount parameter by subjecting the intake 
air amount parameter to predetermined filtering 
processing, and the identification step includes 
identifying the model parameters according to the 
generated filtered value of the intake air amount 
parameter and the plurality of simulation values. 

With the arrangements of these preferred 
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embodiments, it is possible to obtain the same 
advantageous effects as provided by the corresponding 
preferred embodiments of the third aspect of the 
present invention - 

To attain the first object, in a seventh aspect 
of the present invention, there is provided a control 
unit including a control program for causing a computer 
to control a plant, wherein the control program causes 
the computer to detect a detection value reflecting a 
behavior of a first internal variable of the plant, 
generate a simulation value simulating the behavior of 
the first internal variable, estimate an estimation 
value of the detection value based on a model defining 
a relationship between the estimation value and the 
simulation value, identify a model parameter of the 
model according to the detected detection value and the 
generated simulation value, such that the estimated 
estimation value becomes equal to the detected 
detection value, and determine a first input to be 
inputted to the plant, according to the identified 
model parameter. 

With the arrangement of the control unit 
according to the seventh aspect of the present 
invention, it is possible to obtain the same 
advantageous effects as provided by the first aspect of 
the present invention. 

Preferably, the control program causes the 
computer to determine a second input to be inputted to 
the plant such that the detection value is caused to 
converge to a predetermined target value, the first 
internal variable comprising a plurality of first 
internal variables, the simulation value comprising a 
plurality of simulation values simulating respective 
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behaviors of the plurality of first internal variables, 
the model parameter comprising a plurality of model 
parameters ; when the control program causes the 
computer to identify the model parameter, the control 
program causes the computer to identify the plurality 
of model parameters according to the detection value 
and the plurality of simulation values such that the 
estimated estimation value becomes equal to the 
detected detection value; and when the control program 
causes the computer to determine the first input, the 
control program causes the computer to determine the 
first input such that the identified model parameters 
converge to an average value thereof. 

Preferably, when the control program causes the 
computer to determine the first input, the control 
program causes the computer to calculate a learned 
correction value of the first input, using a sequential 
statistical algorithm, correct the first input using 
the calculated learned correction value, and input the 
corrected first input to the plant. 

More preferably, when the control program causes 
the computer to calculate the learned correction value, 
the control program causes the computer to calculate 
the learned correction value of the first input using a 
regression equation in which the learned correction 
value is used as a dependent variable and a second 
internal variable having influence on the first 
internal variable is used as an independent variable, 
and calculate a regression coefficient and a constant 
term of the regression equation with the sequential 
statistical algorithm . 

Preferably, when the control program causes the 
computer to determine the first input, the control 
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program causes the computer to determine an input 
component contained in the first input based on a 
difference between the model parameter and a 
predetermined target value. 

More preferably, when the control program causes 
the computer to determine the first input, the control 
program causes the computer to determine other input 
components than the input component contained in the 
first input, based on the model parameter. 

Preferably, when the control program causes the 
computer to determine the first input, the control 
program causes the computer to determine the first 
input according to the model parameter with a response- 
specified control algorithm. 

Preferably, when the control program causes the 
computer to identify the model parameter, the control 
program causes the computer to identify the model 
parameter by a fixed gain method. 

Further preferably, when the control program 
causes the computer to identify the model parameter, 
the control program causes the computer to identify the 
model parameter by calculating a model parameter 
reference value according to the second internal 
variable, and add a predetermined correction component 
to the calculated model parameter reference value. 

Preferably, the control program causes the 
computer to delay one of the detection value and the 
simulation value by a predetermined delay time period, 
and when the control program causes the computer to 
identify the model parameter, the control program 
causes the computer to identify the model parameter 
according to the delayed one of the detection value and 
the simulation value, and the other of the detection 



50 



H02-2749US 



value and the simulation value. 

Preferably, the control program causes the 
computer to generate a filtered value of the detection 
value by subjecting the detection value to 
predetermined filtering processing, and when the 
control program causes the computer to identify the 
model parameter, the control program causes the 
computer to identify the model parameter according to 
the filtered value of the detection value and the 
simulation value. 

With the arrangements of these preferred 
embodiments, it is possible to obtain the same 
advantageous effects as provided by the corresponding 
preferred embodiments of the first aspect of the 
present invention . 

To attain the second object, in an eighth aspect 
of the present invention, there is provided a control 
unit for an internal combustion engine including a 
plurality of cylinders, a plurality of exhaust passages 
extending from the plurality of cylinders, respectively, 
and one exhaust passage into which the plurality of 
exhaust passages are combined, the control unit 
including a control program for causing a computer to 
perform a control process for controlling an amount of 
fuel to be supplied to the plurality of cylinders, on a 
cylinder-by- cylinder basis, thereby controlling an air- 
fuel ratio of exhaust gases emitted from the plurality 
of cylinders, wherein the control program causes the 
computer to determine an amount of fuel to be supplied 
to each of the plurality of cylinders, detect an air- 
fuel ratio parameter indicative of an air-fuel ratio of 
exhaust gases in the one exhaust passage, generate a 
plurality of simulation values simulating respective 
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behaviors of air-fuel ratios of exhaust gases emitted 
from the plurality of cylinders, estimate an estimation 
value of the air-fuel ratio parameter based on a model 
defining a relationship between the estimation value 
and the plurality of simulation values, identify a 
plurality of model parameters of the model according to 
the detected air-fuel ratio parameter and the generated 
plurality of simulation values, such that the 
estimation value of the air-fuel ratio parameter 
becomes equal to the detected air-fuel ratio parameter, 
calculate a first correction value for correcting the 
amount of fuel to be supplied to the plurality of 
cylinders, according to the identified plurality of 
model parameters, on a cylinder-by-cylinder basis, and 
correct the determined amount of fuel according to the 
calculated first correction value, on a cylinder-by- 
cylinder basis. 

With the arrangement of the control unit 
according to the eighth aspect of the present invention, 
it is possible to obtain the same advantageous effects 
as provided by the second aspect of the present 
invention . 

Preferably, the control program causes the 
computer to calculate a second correction value for 
correcting the amount of fuel to be supplied to each 
cylinder, such that the air- fuel ratio parameter is 
caused to converge to a predetermined target value, and 
correct the amount of fuel to be supplied to each 
cylinder according to the calculated second correction 
value, and when the control program causes the computer 
to calculate the first correction value, the control 
program causes the computer to calculate the first 
correction value, on a cylinder-by-cylinder basis, such 
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that the identified plurality of model parameters 
converge to an average value thereof. 

Preferably, the control program causes the 
computer to calculate a learned correction value of the 
first correction value with a sequential statistical 
algorithm, on a cylinder-by-cylinder basis, and when 
the control program causes the computer to correct the 
amount fuel, the control program causes the computer to 
correct the amount of fuel further according to the 
calculated learned correction value, on a cylinder-by- 
cylinder basis . 

More preferably, the control program causes the 
computer to detect an operating condition parameter 
indicative of an operating condition of the engine, and 
when the control program causes the computer to 
calculate the learned correction value, the control 
program causes the computer to calculate the learned 
correction value using a regression equation in which 
the learned correction value is used as a dependent 
variable and the detected operating condition parameter 
is used as an independent variable, and calculate a 
regression coefficient and a constant term of the 
regression equation with the sequential statistical 
algorithm. 

Preferably, when the control program causes the 
computer to calculate the first correction value, the 
control program causes the computer to calculate a 
correction value component contained in the first 
correction value based on a difference between the 
identified model parameters and a predetermined target 
value . 

More preferably, when the control program causes 
the computer to calculate the first correction value. 
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the control program causes the computer to calculate 
other correction value components than the correction 
value component contained in the first correction value, 
based on the identified model parameters. 

Preferably, when the control program causes the 
computer to calculate the first correction value, the 
control program causes the computer to calculate the 
first correction value according to the model 
parameters with a response-specified control algorithm. 

Preferably, when the control program causes the 
computer to identify the model parameters of the model, 
the control program causes the computer to identify the 
model parameters by a fixed gain method. 

Further preferably, when the control program 
causes the computer to identify the model parameters of 
the model, the control program causes the computer to 
identify the model parameters by calculating respective 
model parameter reference values according to the 
operating condition parameter, and adding predetermined 
correction components to the calculated model parameter 
reference values, respectively. 

Preferably, the control program causes the 
computer to delay the air- fuel ratio parameter by a 
predetermined delay time period, and when the control 
program causes the computer to identify the model 
parameters of the model, the control program causes the 
computer to identify the model parameters according to 
the delayed air-fuel ratio parameter and the plurality 
of simulation values . 

With the arrangements of these preferred 
embodiments, it is possible to obtain the same 
advantageous effects as provided by the corresponding 
preferred embodiments of the second aspect of the 
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present invention . 

To attain the second object, in a ninth aspect of 
the present invention, there is provided a control unit 
for an internal combustion engine including one intake 
passage, a plurality of intake passages branching from 
the one intake passage, and a plurality of cylinders 
connected to the plurality of intake passages extend, 
respectively, the control unit including a control 
program for causing a computer to perform a control 
process for controlling an amount of fuel to be 
supplied to the plurality of cylinders, on a cylinder- 
by-cylinder basis, thereby controlling an air-fuel 
ratio of exhaust gases emitted from the plurality of 
cylinders, wherein the control program causes the 
computer to determine an amount of fuel to be supplied 
to each of the plurality of cylinders, detect an intake 
air amount parameter indicative of an amount of intake 
air in the one intake passage, generate a plurality of 
simulation values simulating respective behaviors of 
amounts of intake air to be drawn into the plurality of 
cylinders, estimate an estimation value of the intake 
air amount parameter based on a model defining a 
relationship between the estimation value and the 
plurality of simulation values, identifying a plurality 
of model parameters of the model according to the 
detected intake air amount parameter and the generated 
plurality of simulation values, such that the 
estimation value of the intake air amount parameter 
becomes equal to the detected intake air amount 
parameter, calculate a third correction value for 
correcting the amount of fuel to be supplied to the 
plurality of cylinders, according to the identified 
plurality of model parameters, on a cylinder-by- 
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cylinder basis, and correct the determined amount of 
fuel according to the calculated third correction value, 
on a cylinder-by-cylinder basis. 

With the arrangement of the control unit 
according to the ninth aspect of the present invention, 
it is possible to obtain the same advantageous effects 
as provided by the third aspect of the present 
invention . 

Preferably, the engine includes a plurality of 
exhaust passages extending from the plurality of 
cylinders, respectively, and one exhaust passage into 
which the plurality of exhaust passages are combined; 
the control program causes the computer to detect an 
air-fuel ratio parameter indicative of an air-fuel 
ratio of exhaust gases in the one exhaust passage, 
calculate a fourth correction value for correcting the 
amount of fuel to be supplied to each cylinder, such 
that the detected air-fuel ratio parameter is caused to 
converge to a predetermined target value, and correct 
the amount of fuel to be supplied to each cylinder 
according to the calculated fourth correction value; 
and when the control program causes the computer to 
calculate the third correction value, the control 
program causes the computer to calculate the third 
correction value, on a cylinder-by-cylinder basis, such 
that the identified plurality of model parameters 
converge to an average value thereof . 

Preferably, the control program causes the 
computer to calculate a learned correction value of the 
third correction value with a sequential statistical 
algorithm, on a cylinder-by-cylinder basis, and when 
the control program causes the computer to correct the 
amount of fuel, the control program causes the computer 
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to correct the amount of fuel further according to the 
calculated learned correction value, on a cylinder-by- 
cylinder basis. 

More preferably, the control program causes the 
computer to detect an operating condition parameter 
indicative of an operating condition of the engine, and 
when the control program causes the computer to 
calculate the learned correction value, the control 
program causes the computer to calculate the learned 
correction value using a regression equation in which 
the learned correction value is used as a dependent 
variable and the detected operating condition parameter 
is used as an independent variable, and calculate a 
regression coefficient and a constant term of the 
regression equation with the sequential statistical 
algorithm. 

Preferably, when the control program causes the 
computer to calculate the third correction value, the 
control program causes the computer to calculate a 
correction value component contained in the third 
correction value based on a difference between the 
identified model parameters and a predetermined target 
value . 

More preferably, when the control program causes 
the computer to calculate the third correction value, 
the control program causes the computer to calculate 
other correction value components than the correction 
value component contained in the third correction value, 
based on the identified model parameters. 

Preferably, when the control program causes the 
computer to calculate the third correction value, the 
control program causes the computer to calculate the 
third correction value according to the model 
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parameters with a response-specified control algorithm. 

Preferably, when the control program causes the 
computer to identify the model parameters of the model, 
the control program causes the computer to identify the 
model parameters by a fixed gain method. 

Further preferably, when the control program 
causes the computer to identify the model parameters of 
the model, the control program causes the computer to 
identify the model parameters by calculating respective 
model parameter reference values according to the 
operating condition parameter, and adding predetermined 
correction components to the calculated model parameter 
reference values, respectively. 

Preferably, the control program causes the 
computer to delay the plurality of simulation values by 
a predetermined delay time period, and when the control 
program causes the computer to identify the model 
parameters of the model, the control program causes the 
computer to identify the model parameters according to 
the intake air amount parameter and the delayed 
plurality of simulation values. 

Preferably, the control program causes the 
computer to generate a filtered value of the intake air 
amount parameter by subjecting the intake air amount 
parameter to predetermined filtering processing, and 
when the control program causes the computer to 
identify the model parameters of the model, the control 
program causes the computer to identify the model 
parameters according to the generated filtered value of 
the intake air amount parameter and the plurality of 
simulation values. 

With the arrangements of these preferred 
embodiments, it is possible to obtain the same 
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advantageous effects as provided by the corresponding 
preferred embodiments of the third aspect of the 
present invention . 

The above and other objects, features, and 
advantages of the invention will become more apparent 
from the following detailed description taken in 
conjunction with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram schematically showing 
the arrangement of a control system according to a 
first embodiment of the present invention and an 
internal combustion engine to which the control system 
is applied; 

FIG. 2 is a block diagram showing the arrangement 
of the control system according to the first 
embodiment ; 

FIG. 3 is a schematic diagram useful in 
explaining a calculation algorithm with which an air- 
fuel ration variation coefficient <E> is calculated by 
an adaptive observer of a first air-fuel ratio 
controller; 

FIG. 4 is a diagram showing mathematical 
expressions of the calculation algorithm with which the 
air-fuel ration variation coefficient O is calculated 
by the adaptive observer; 

FIG. 5 is a block diagram showing the 
configuration of the adaptive observer; 

FIG. 6 is a diagram showing changes in a 
simulation value KACT_OS output from a signal generator 
of the adaptive observer; 

FIG. 7 is a diagram showing mathematical 
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expressions of an I-PD control algorithm with which an 
air-fuel ratio variation correction coefficient KOBSV 
is calculated; 

FIG. 8A is a diagram showing the relationship 
between an exhaust gas volume ESV and the air- fuel 
ratio variation correction coefficient KOBSV; 

FIG. 8B is a diagram showing the relationship 
between the exhaust gas volume ESV, the air-fuel ratio 
variation correction coefficient KOBSV, and a learned 
correction value KOBSV_LS thereof; 

FIG. 9 is a diagram showing mathematical 
expressions of a calculation algorithm with which the 
learned correction value KOBSV_LSi of the air- fuel 
ratio variation correction coefficient is calculated; 

FIG. 10 is a diagram showing mathematical 
expressions useful for explaining a calculation 
algorithm with which a feedback correction coefficient 
KSTR is calculated by a second air-fuel ratio 
controller; 

FIG. 11 is a diagram showing mathematical 
expressions of a calculation algorithm with which the 
feedback correction coefficient KSTR is calculated by 
the second air-fuel ratio controller; 

FIG. 12 is a flowchart showing an air-fuel ratio 
control process; 

FIG. 13 is a flowchart showing a process for 
calculating a model parameter vector 9 , executed in a 
step 6 in FIG. 12; 

FIG. 14 is a flowchart showing a KSTR-calculating 
process executed in a step 7 in FIG. 12; 

FIG. 15 is a flowchart showing a process executed 
in a step 8 in FIG. 12, for calculating a vector (f> of 
the air-fuel ratio variation coefficient; 
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FIG. 16 is a flowchart showing a process for 
calculating the air-fuel ratio variation correction 
coefficient KOBSVj., executed in a step 9 in FIG. 12; 

FIG. 17 is a flowchart showing a process for 
calculating the learned correction value KOBSV_LSi of 
the air-fuel ratio variation correction coefficient, 
executed in a step 10 in FIG. 12; 

FIG. 18 is a timing chart showing an example of 
operations of the air-fuel ratio control executed by 
the control system according to the first embodiment; 

FIG. 19 is a timing chart showing a comparative 
example of operations of the air-fuel ratio control; 

FIG. 20 is a diagram showing respective groups of 
mathematical expressions of a PID control algorithm, an 
IP-D control algorithm, and a response- specif led 
control algorithm with each of which the air-fuel ratio 
variation correction coefficient KOBSVi is calculated; 

FIG. 21 is a diagram showing an identification 
algorithm with which the vector (£> of the air- fuel 
ratio variation coefficient is identified by a fixed 
gain method to which a 6 correcting method is applied; 
/ FIG. 2 2 is a diagram showing an example of a 

table used in calculating reference values Obasei 

FIG. 2 3 is a block diagram showing the 
arrangement of components of a control system according 
to a second embodiment of the present invention, with a 
third air-fuel ratio controller as the center thereof; 

FIG. 24 is a block diagram showing the 
arrangement of components of the control system 
according to the second embodiment, with first and 
second air-fuel ratio controllers as the center 
thereof ; 

FIG. 25 is a diagram showing a waveform of 
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pulsation of intake air detected by an air flow sensor; 

FIG. 26 is a schematic diagram useful for 
explaining a calculation algorithm with which an intake 
air amount variation coefficient V is calculated by an 
adaptive observer of the third air-fuel ratio 
controller; 

FIG. 27 is a diagram showing mathematical 
expressions of the calculation algorithm with which the 
intake air amount variation coefficient V is 
calculated by the adaptive observer; 

FIG. 28 is a block diagram showing the 
configuration of the adaptive observer; 

FIG. 29 is a diagram showing changes in a 
simulation value GAIR_OS output from a signal generator 
of the adaptive observer; 

FIG. 30 is a diagram showing mathematical 
expressions of an I-PD control algorithm with which an 
intake air amount variation correction coefficient 
KICYL is calculated; 

FIG. 31 is a diagram showing the relationship 
between the exhaust gas volume ESV, the intake air 
amount variation correction coefficient KICYL , and a 
learned correction value KICYL_LS thereof; 

FIG. 32 is a diagram showing mathematical 
expressions of a calculation algorithm with which the 
learned correction value KICYL_LS of the intake air 
amount variation correction coefficient is calculated; 

FIG. 33 is a flowchart showing an air-fuel ratio 
control process executed by the control system 
according to the second embodiment; 

FIG. 34 is a flowchart showing a process for 
calculating a vector i/> of the intake air amount 
variation coefficient, executed in a step 111 in FIG. 
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34; 

FIG. 35 is a flowchart showing a process for 
calculating an intake air amount variation correction 
coefficient KICYL i# executed in a step 112 in FIG. 34; 

FIG. 36 is a flowchart showing a process for 
calculating the learned correction value KICYL_LS ± of 
the intake air amount variation correction coefficient, 
executed in a step 113 in FIG. 34; 

FIG. 37 is a block diagram showing the 
configuration of a variation of the adaptive observer 
of the third air-fuel ratio controller; 

FIG. 38 is a diagram useful in explaining an 
example of a filter of the variation of the adaptive 
observer; 

FIG. 39 is a diagram useful in explaining another 
example of the filter of the variation of the adaptive 
observer; 

FIG. 40 is a diagram showing mathematical 
expressions of the filter of the variation of the 
adaptive observer, and mathematical expressions of a 
calculation algorithm of the intake air amount 
variation coefficient Y ; 

FIG. 41 is a diagram showing respective groups of 
mathematical expressions of an IP-D control algorithm 
and a response-specified control algorithm with each of 
which the intake air amount variation correction 
coefficient KICYL ± is calculated; 

FIG. 42 is a diagram showing mathematical 
expressions of an identification algorithm with which 
the vector il> of the intake air amount variation 
coefficient is identified by a fixed gain method to 
which the 6 correction method is applied; and 

FIG. 43 is a diagram showing a table used in the 



H02-2749US 

63 

calculation of reference values ?basei. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

The invention will now be described in detail 
with reference to drawings showing preferred 
embodiments thereof. Referring first to FIG. 1, there 
is schematically shown the arrangement of an control 
system 1 according to a first embodiment of the present 
invention and an internal combustion engine 3 , as a 
plant, to which the control system 1 is applied. As 
shown in FIG. 1, the control system 1 includes an ECU 2 
which controls the amount of fuel injected into the 
internal combustion engine (hereinafter simply referred 
to as "the engine") 3 according to operating conditions 
of the engine 3, to thereby control the air-fuel ratio 
of the mixture, as described in detail hereinafter. 

The engine 3 is an inline four-cylinder gasoline 
engine installed on an automotive vehicle, not shown, 
and has first to fourth cylinders #1 to #4 (a plurality 
of cylinders). The engine 3 has an intake pipe 4 which 
includes a main pipe 4a (one intake air passage), and 
an intake manifold 4b connected thereto. A throttle 
valve 5 is arranged across an intermediate portion of 
the main pipe 4a. 

. At respective locations upstream and downstream 
of the throttle valve 5, there are arranged an air flow 
sensor 9 and an intake pipe absolute pressure sensor 11. 
The air flow sensor 9 detects the amount of intake air 
GAIR (detection value, intake air amount parameter) 
drawn into the engine via the intake pipe 4 , and 
delivers a signal indicative of the detected amount of 
intake air to the ECU 2. 
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Further, the intake pipe absolute pressure sensor 
11 is implemented e.g. by a semiconductor pressure 
sensor, which detects the intake pipe absolute pressure 
PBA (detection value, intake air amount parameter) of 
the intake pipe 4, and delivers a signal indicative of 
the detected intake pipe absolute pressure PBA to the 
ECU 2. In the present embodiment, the air flow sensor 
9 forms the detection means, the operating condition 
parameter-detecting means, and the intake air amount 
parameter-detecting means, while the intake pipe 
absolute pressure sensor 11 forms the detection means 
and the intake air amount parameter-detecting means. 

In the vicinity of the throttle valve 5 disposed 
in the main pipe 4a, there is provided a throttle valve 
opening sensor 10 implemented e.g. by a potentiometer, 
for detecting the degree of opening (hereinafter 
referred to as "throttle valve opening") TH of the 
throttle valve 5 and delivering an electric signal 
indicative of the sensed throttle valve opening TH to 
the ECU 2. 

The intake manifold 4b of the intake pipe 4 is 
comprised of a collecting section 4c (one intake 
passage) connected to the main pipe 4a, and four branch 
portions 4d (plurality of intake passages) branching 
from the collecting section 4c and connected to the 
four cylinders #1 to #4, respectively. In the branch 
portions 4d, injectors 6 are inserted at respective 
locations upstream of intake ports, not shown, for the 
cylinders. During operation of the engine 3, each 
injector 6 is controlled in respect of a fuel injection 
amount, i.e. a time period over which the injector 6 is 
open, and fuel injection timing, by a drive signal 
delivered from the ECU 2. 
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Further, an engine coolant temperature sensor 12 
implemented e.g. by a thermistor is mounted in the 
cylinder block of the engine 3. The engine coolant 
temperature sensor 12 senses an engine coolant 
temperature TW which is the temperature of an engine 
coolant circulating through the cylinder block of the 
engine 3 and delivering a signal indicative of the 
sensed engine coolant temperature TW to the ECU 2 . 

A crank angle position sensor 13 (operating 
condition parameter-detecting means) is provided for a 
crankshaft, not shown, of the engine 3, for delivering 
a CRK signal and a TDC signal, which are both pulse 
signals, to the ECU 2 in accordance with rotation of 
the crankshaft. 

Each pulse of the CRK signal is generated 
whenever the crankshaft rotates through a predetermined 
angle (e.g. 30 degrees). The ECU 2 determines a 
rotational speed (hereinafter referred to as "the 
engine speed") NE of the engine 3, based on the CRK 
signal. The TDC signal indicates that each piston, not 
shown, in an associated cylinder is in a predetermined 
crank angle position immediately before the TDC 
position at the start of the intake stroke, and each 
pulse of the TDC signal is generated whenever the 
crankshaft rotates through a predetermined angle. 

On the other hand, the exhaust pipe 7 includes an 
exhaust manifold 7b connected to the four cylinders #1 
to #4, and an main pipe 7a connected to a collecting 
section 7j of the exhaust manifold 7b. The exhaust 
manifold 7b is configured such that four exhaust pipe 
sections 7c to 7f (a plurality of exhaust passages) 
extending from the four cylinders #1 to #4 are combined 
into two collecting sections, and the two collecting 
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sections are combined into one collecting section. 
That is, the exhaust manifold 7b is comprised of two 
exhaust pipe sections 7c and 7f extending from the 
respective first and fourth cylinders #1 and #4, a 
collecting section 7g into which these exhaust pipe 
sections 7c and 7f are combined, two exhaust pipe 
sections 7d and 7e extending from the respective second 
and third cylinders #2 and #3, and a collecting section 
7h into which these exhaust pipe sections 7d and 7e are 
combined, and a collecting section 7j (one exhaust 
passage) into which the two collecting sections 7g and 
7h are combined, all of these components being 
integrally formed in one piece. Due to such a 
configuration, the exhaust manifold 7b has a lower 
resistance to the flow of exhaust gases than a 
conventional exhaust manifold in which four exhaust 
pipe sections are directly combined into one collecting 
section. This enables the engine 3 to deliver higher 
power output and higher torque, compared with those 
having the conventional exhaust manifold. 

A first catalytic device 8a and a second 
catalytic devices 8b are arranged in the exhaust pipe 7 
from upstream to downstream in the mentioned order in a 
spaced relationship at respective locations of the main 
pipe 7a of the intake pipe 7. Each catalytic device 8 
is a combination of a NOx catalyst and a three-way 
catalyst, and the NOx catalyst is comprised of a 
honeycomb structure base, an iridium catalyst (sintered 
body of silicon carbide whiskers carrying iridium and 
silica) coated on the surface of the honeycomb 
structure base, and Perovskite double oxide (sintered 
body of LaCo0 3 powder and silica) further coated on the 
iridium catalyst. The catalytic device 8 eliminates 
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NOx from exhaust gases emitted during a lean burn 
operation of the engine 3 by oxidation-reduction 
reaction-catalyzing action of the NOx catalyst, and 
eliminates CO, HC, and NOx from exhaust gases emitted 
during other operations of the engine 3 than the lean 
burn operation by oxidation-reduction reaction- 
catalyzing action of the three-way catalyst. 

An oxygen concentration sensor (hereinafter 
referred to as "the 02 sensor") 15 is inserted into the 
main pipe 7a between the first and second catalytic 
devices 8a and 8b. The 02 sensor 15 is comprised of a 
zirconia layer and platinum electrodes, and delivers to 
the ECU 2 an output Vout dependent on the concentration 
of oxygen contained in exhaust gases downstream of the 
first catalytic device 8a. The output Vout assumes a 
high-level voltage value (e.g. 0.8 V) when an air- fuel 
mixture having a richer air- fuel ratio than the 
stoichiometric air-fuel ratio has been burned, whereas 
it assumes a low- level voltage value (e.g. 0.2V) when 
an air-fuel mixture having a leaner air-fuel ratio than 
the stoichiometric air-fuel ratio has been burned. 
Further, when the air- fuel ratio of the mixture is 
close to the stoichiometric air-fuel ratio, the output 
Vout assumes a predetermined target value Vop (e.g. 0.6 
V) between the high-level and low-level voltage values. 

Further, a LAF sensor 14 is mounted in the 
vicinity of the collecting section 7d of the exhaust 
manifold 7a. The LAF sensor 14 (detection means, air- 
fuel ratio parameter-detecting means) is formed by 
combining a sensor similar to the 02 sensor 15 and a 
detection circuit, such as a linearizer, and detects 
the concentration of oxygen contained in exhaust gases 
linearly over a wide range of the air-fuel ratio 
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ranging from a rich region to a lean region, thereby 
delivering an output proportional to the sensed oxygen 
concentration to the ECU 2. The ECU 2 calculates a 
detected air-fuel ratio KACT (detection value, air-fuel 
ratio parameter) indicative of the air-fuel ratio of 
exhaust gases at the collecting section 7j based on the 
output from the LAF sensor 14. The detected air-fuel 
ratio KACT is expressed as an equivalent ratio 
proportional to the reciprocal of the air-fuel ratio. 

Further, the ECU 2 has an accelerator pedal 
opening sensor 16; an atmospheric pressure sensor 17, 
an intake air temperature sensor 18, and a vehicle 
speed sensor 19, connected thereto. The accelerator 
pedal opening sensor 16 detects a depression amount 
(hereinafter referred to as "the accelerator pedal 
opening") AP of an accelerator pedal, not shown, of the 
vehicle and delivers a signal indicative of the sensed 
accelerator pedal opening AP to the ECU 12. Further, 
the atmospheric pressure sensor 17, the intake air 
temperature sensor 18, and the vehicle speed sensor 19 
detect atmospheric pressure AP, intake air temperature 
TA, and vehicle speed VP, respectively, and delivers 
respective signals indicative of the detected 
atmospheric pressure AP, intake air temperature TA, and 
vehicle speed VP to the ECU 2 . 

Next, the ECU 2 will be described. The ECU 2 is 
implemented by a microcomputer including an 
input /output interface, a CPU, a RAM, and a ROM, none 
of which is shown. The ECU 2 determines operating 
conditions of the engine 3, based on the outputs from 
the aforementioned sensors 9 to 19. Further, the ECU 2 
executes an air-fuel ratio control process, which will 
be described in detail hereinafter, according to 
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control programs read from the ROM, using data stored 
in the RAM, and the like, to thereby calculate a target 
air-fuel ratio KCMD, a feedback correction coefficient 
KSTR, an air-fuel ratio variation correction 
coefficient KOBSV ± , and a learned correction value 
KOBSV_LSi thereof. Furthermore, as described 
hereinafter, based on these calculated values of KCMD, 
KSTR, KOBSVi, KOBSV_LS , , and so forth, the ECU 2 
calculates a final fuel injection amount TOUT ± for each 
injector 6, on a cylinder-by-cylinder basis, and drives 
the injector 6 by a drive signal generated based on the 
calculated final fuel injection amount TOUTi , to 
thereby control the air-fuel ratio of the mixture, i.e. 
air-fuel ratio of exhaust gases, on a cylinder-by- 
cylinder basis. It should be noted that the subscript 
"t" in TOUTi represents a cylinder number indicative of 
a number assigned to each cylinder (i = 1 to 4), and 
this also applies to the aforementioned air-fuel ratio 
variation correction coefficient KOBSVi, and the 
learned correction value KOBSV_LSs, and parameters 
referred to hereinafter. 

It should be noted that in the present embodiment , 
the ECU 2 forms simulation value-generating means, 
estimation means, identification means, first control 
means, second control means, learned correction value- 
calculating means, correction means, input means, delay 
means, fuel amount -determining means, first correction 
value-calculating means, first fuel amount-correcting 
means, second correction value-calculating means, 
second fuel amount -correcting means, operating 
condition parameter-detecting means, third correction 
value-calculating means, third fuel amount -correcting 
means, fourth correction value-calculating means, and 
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fourth fuel amount -correcting means. 

As shown in FIG. 2, the control system 1 is 
comprised of a basic fuel injection amount-calculating 
section 20 , a first air-fuel ratio controller 30, a 
second air-fuel ratio controller 40, and a fuel 
attachment -dependent correcting section 50, which are 
all implemented by the ECU 2. In the control system 1, 
the basic fuel injection amount -calculating section 20, 
as the fuel amount -determining means, calculates a 
basic fuel injection amount TIBS according to the 
intake air amount GAIR by searching a table, not shown. 

Further, as described in detail hereinafter, to 
correct variation in air-fuel ratio between the 
cylinders, the first air-fuel ratio controller 30 
calculates the air-fuel ratio variation correction 
coefficient KOBSV ± and the learned correction value 
KOBSV_LSi thereof, and the second air-fuel ratio 
controller 40 calculates the feedback correction 
coefficient KSTR so as to cause the detected air-fuel 
ratio KACT to converge to the target air- fuel ratio 
KCMD. Then, a demanded fuel injection amount TCYL ± is 
calculated on a cylinder-by-cylinder basis by 
multiplying the basic fuel injection amount TIBS by a 
corrected target air-fuel ratio KCMDM, a total 
correction coefficient KTOTAL, the feedback correction 
coefficient KSTR, the air-fuel ratio variation 
correction coefficient KOBSV i# and the learned 
correction value KOBSV_LS ± thereof . Then, the fuel 
attachment -dependent correcting section 50 calculates 
the final fuel injection amount TOUT ± based on the 
demanded fuel injection amount TCYL ± , on a cylinder-by- 
cylinder basis. 

Next, a description will be given of the first 
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air-fuel ratio controller 30. The first air-fuel ratio 
controller 30 (first fuel amount -correcting means )is 
for correcting variation in air-fuel ratio between the 
cylinders, and is comprised of an adaptive observer 31, 
an air-fuel ratio variation correction coefficient- 
calculating section 32, a learned correction value- 
calculating section 33, and a multiplication section 
34. 

In this first air-fuel ratio controller 30, with 
algorithms, described hereinafter, the adaptive 
observer 31 (simulation value-generating means, 
estimation means, identification means, delay means) 
calculates an air-fuel ratio variation coefficient 0 ±/ 
on a cylinder-by-cylinder basis, and the air-fuel ratio 
variation correction coefficient-calculating section 32 
(first control means, first correction value- 
calculating means) calculates the air-fuel ratio 
variation correction coefficient KOBSV ± , on a cylinder- 
by-cylinder basis. Further, the learned correction 
value-calculating section 33 (learned correction value- 
calculating means) calculates the learned correction 
value KOBSV_LSi of the air- fuel ratio variation 
correction coefficient, on a cylinder-by-cylinder basis. 
Further, the multiplication section 34 (correction 
means) multiplies the air-fuel ratio variation 
correction coefficient KOBSVi to KOBSV 4 by the learned 
correction value K0BSV_LSi to KOBSV_LS 4 , respectively. 
That is, the air-fuel ratio variation correction 
coefficient KOBSV ± is corrected by the learned 
correction value KOBSV_LS ± . 

Next, a description will be given of an algorithm 
of the adaptive observer 31. First, as shown in FIG. 3, 
the intake system of the engine 3 is regarded as a 
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system which is represented by four simulation values 
KACT_OSi to KACT_OS 4 and four air- fuel ratio variation 
coefficients $1 to <3> 4 . These simulation values 
KACT_OS ± are values simulating the exhaust timing and 
exhaust behavior of exhaust gases, on a cylinder-by- 
cylinder basis, and the air-fuel ratio variation 
coefficient <E>i represents variation in air-fuel ratio 
of exhaust gases between the cylinder and the amount of 
change in the exhaust behavior. When this system is 
modeled into a discrete time system model, an equation 
(1) shown in FIG. 4 can be obtained. In the equation 
(1), the symbol k represents a discretized time, and 
indicates that each discrete data (time-series data) 
with (k) is data sampled whenever a pulse of the TDC 
signal is generated. This also applies to discrete 
data referred to hereinafter. (Each discrete data may 
be data sampled whenever a pulse of the CRK signal is 
generated.) Further, the symbol d represents dead time 
(predetermined delay time) which the exhaust gases 
emitted from each cylinder take to reach the LAF sensor 
LAF, and is set to a predetermined fixed value in the 
present embodiment in advance. The dead time d may be 
set depending on an operating condition (engine speed 
NE) of the engine 3. 

The adaptive observer 31 according to the present 
embodiment uses an equation formed by replacing the 
left side of the equation (1) by the estimation value 
KACT_EST of the air-fuel ratio, i.e. an equation (2) in 
FIG. 4 as a model, and a signal generator 31a generates 
the simulation value KACT_OSi. At the same time, the 
vector <b (k) of the air-fuel ratio variation 
coefficient 0 ± as a model parameter of the equation 



73 



H02-2749US 



(2) is identified by a variable-gain sequential least- 
squares method expressed by equations (3) to (9) in FIG. 
4 such that the estimation value KACT_EST(k) becomes 
equal to the detected air-fuel ratio KACT(k) . 

The symbol KP(k) in the equation (3) represents a 
vector of a gain coefficient, and the symbol ide(k) 
represents an identification error. Further, (f>(k) T in 
the equation (4) represents a transposed matrix of <t> 
(k). It should be noted in the following description, 
the notation of "vector" is omitted unless otherwise 
required. The identification error ide(k) in the 
equation (3) is calculated using the equations (5) to 
(7) in FIG. 4, and the symbol £ (k) in the equation (6) 
represents a vector of the simulation values defined by 
the equation (7). Further, the vector KP(k) of the 
gain coefficient is calculated by an equation (8) in 
FIG. 4, and the symbol P(k) in this equation (8) 
represents a square matrix of order 4 defined by an 
equation (9) in FIG. 4. 

This adaptive observer 31 identifies the vector 
^ (k) of the air-fuel ratio variation coefficient < $>± 
with the algorithm based on the sequential least - 
squares method shown in the equations (2) to (9). This 
makes it possible to remove (filter off) noise-like 
fluctuating components of the exhaust behavior caused 
by a sudden change in the operating condition of the 
engine 3 from the air-fuel ratio variation coefficient 
$i, and thereby calculate the air-fuel ratio variation 
coefficient <& ± as a value substantially indicative of 
variation in air-fuel ratio between the cylinders. 

The configuration of the adaptive observer 31 can 
be represented by a block diagram shown in FIG. 5. 
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That is, as shown in FIG. 5, in the adaptive observer 
31 , the signal generator 31a generates the vector £(k) 
of the simulation values KACT_OSi. More specifically, 
as shown in FIG- 6, the signal generator 31a generates 
the simulation values KACT_OS ± as signal values each 
having a waveform of a combination of triangular waves 
and trapezoidal waves formed such that the total sum of 
the simulation values constantly becomes equal to a 
value of 1. Further, the multiplier 31b generates the 
estimation value KACT__EST(k) of the air-fuel ratio as a 
value obtained by multiplying the vector £ (k) of the 
simulation values by the vector (f>(k-l) of the air-fuel 
ratio variation coefficient. Then, the differentiator 
31d generate the identification error ide(k) as the 
difference between the detected air-fuel ratio KACT(k) 
and the estimation value KACT_EST(k) . 

Further, a logic unit 31e generates the vector 
KP(k) of the gain coefficient based on the vector £ (k) 
of the simulation values, and a multiplier 31f 
generates the product [ide(k) • KP(k) ] of the 
identification error ide(k) and the vector KP(k) of the 
gain coefficient. Next, an adder 31g generates the 
vector <t> (k) of the air-fuel ratio variation 
coefficient as the sum of the product [ide(k) -KP(k)] 

and the delayed vector (f>(k-l) of the air-fuel ratio 
variation coefficient . 

Next, a description will be given of an algorithm 
with which the air-fuel ratio variation correction 
coefficient-calculating section 32 calculates the air- 
fuel ratio variation correction coefficient 
KOBSVi(first input, first correction value). In the 
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air-fuel ratio variation correction coefficient - 
calculating section 32, first, the moving average value 
( E>ave(k) of the air-fuel ratio coefficient is 
calculated based on the air-fuel ratio variation 
coefficient 0 ± (k) calculated by the adaptive observer 
31, on a cylinder-by-cylinder basis, by an equation 
(10) in FIG. 7. Next, the air-fuel ratio variation 
correction coefficient KOBSVi is calculated with I-PD 
control (proportional/differential-preceding PID 
control) algorithm, on a cylinder-by-cylinder basis, 
such that the air-fuel ratio variation coefficient O 
i(k) converges to the moving average value $ave(k) 
thereof. This I-PD control algorithm is expressed by 
equations (11) and (12) in FIG. 7. The symbol e(k) in 
the equation (12) represents a following error. 

As described above, the air-fuel ratio variation 
correction coefficient-calculating section 32 
calculates the air-fuel ratio variation correction 
coefficient KOBSVi, with the I-PD control algorithm, 
such that the air-fuel ratio variation coefficient O 
i(k) converges to the moving average value Oave(k) 
thereof. This is for providing control such that the 
converging behavior of the air- fuel ratio variation 
coefficient ( E > i (k) to the moving average value Oave(k) 
thereof does not suffer from overshooting, thereby 
preventing the air-fuel ratio control by the first air- 
fuel ratio controller 30 for correcting variation in 
air- fuel ratio between the cylinders and the air- fuel 
ratio control by the second air-fuel ratio controller 
40, described hereinafter, from interfering with each 
other . 

Next, a description will be given of an algorithm 
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with which the learned correction value-calculating 
section 33 calculates the learned correction value 
KOBSV_LSi of the air-fuel ratio variation correction 
coefficient KOBSVi. The air- fuel ratio variation 
correction coefficient KOBSVi is susceptible to 
operating conditions of the engine 3, and when the 
operating conditions of the engine 3 are changed, the 
coefficient KOBSVi is changed accordingly. FIG. 8A 
shows the relationship between an exhaust gas volume 
ESV(k) as an operating condition parameter indicative 
of an operating condition of the engine and the air- 
fuel ratio variation correction coefficient KOBSV ± (k) . 
This exhaust gas volume ESV(k) (second internal 
variable, operating condition parameter) is an 
estimation value of the space velocity, and calculated 
using an equation (13) in FIG. 9. It should be noted 
that in the equation (13), the symbol SVPRA represents 
a predetermined coefficient determined in advance by 
the displacement of the engine 3 . 

Referring to FIG. 8A, it can be seen that an 
approximate value i.e. estimation value of the air-fuel 
ratio variation correction coefficient KOBSV t (k) can be 
calculated by a first-degree equation using the air- 
fuel ratio variation correction coefficient KOBSV ± (k) 
as a dependent variable and the exhaust gas volume 
ESV(k) as an independent variable (see FIG. 8B) . 
Therefore, in the learned correction value-calculating 
section 33, the learned correction value KOBSV__LS ± (k) 
of the air-fuel ratio variation correction coefficient 
is defined as an estimation value calculated by a 
regression equation, i.e. an equation (14) in FIG. 9, 
and a vector #OBSV_LSi(k) of a regression coefficient 
AOBSV_LS ± and a constant term BOBSV_LSi (hereinafter 
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referred to as "the regression coefficient vector") is 
calculated by a sequential least -squares method 
expressed by equations (15) to (21) in FIG. 9. 

In this equation (15), the symbol KQi(k) 
represents a vector of a gain coefficient, and the 
symbol EoVi(k) represents an error. Further, the error 
Eovi(k) is calculated using equations (17) to (19) in 
FIG. 9. Further, the vector KQi(k) of the gain 
coefficient is calculated using an equation (20) in FIG. 
9, and the symbol Q±(k) in this equation (20) 
represents a square matrix of order 2 defined by an 
equation (21) in FIG. 9. 

Further, the learned correction value 
KOBSV_LSi(k) is more specifically calculated using an 
equation (22) in FIG. 9. It should be noted that when 
the engine 3 is in an extreme operating condition or 
operating environment, the calculation of the 
regression coefficient AOBSV_LS ± and the constant term 
BOBSV_LSi by the sequential least - squares method is 
avoided, and the immediately preceding value 9 
OBSV_LSi(k-l ) of the regression coefficient vector is 
used as the current value #0BSV_LSi(k) in calculation 
of the learned correction value KOBSV__LS ± (k) . 

With the algorithm expressed by the equations 
(13) and (15) to (22), the learned correction value- 
calculating section 33 calculates the learned 
correction value KOBSV_LS ± (k) such that the learned 
correction value KOBSV_LSi(k) converges to the product 
of the learned correction value KOBSV_LS ± (k) and the 
air-fuel ratio variation correction coefficient 
KOBSVi(k) . 

Next, a description will be given of the second 
air-fuel ratio controller 40. The second air-fuel 
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ratio controller 40 is formed as an STR (Self Tuning 
Regulator) comprised of an onboard identifier 41 and an 
STR controller 42. The second air-fuel ratio 
controller 40 calculates the feedback correction 
coefficient KSTR such that the detected air- fuel ratio 
KACT converges to the target air-fuel ratio KCMD 
(predetermined target value). More specifically, with 
an algorithm described hereinbelow, the onboard 
identifier 41 identifies a model parameter vector 0 of 
the first cylinder #1, and the STR controller 42 
calculates the feedback correction coefficient KSTR 
(second input, second and fourth correction values). 
In the present embodiment, this second air-fuel ratio 
controller 40 forms the second control means, the 
second correction value-calculating means, the second 
fuel amount -correcting means, the fourth correction 
value-calculating means, and the fourth fuel amount- 
correcting means. 

First, the first cylinder #1 is regarded as a 
controlled object to which is inputted the feedback 
correction coefficient KSTR and from which is outputted 
the detected air-fuel ratio KACT, and this controlled 
object is modeled into a discrete- time system model, 
which is expressed by an equation (23) in FIG. 10. In 
the equation (23), the symbol n represents a 
discretized time, and indicates that each discrete data 
with (n) is sampled every combustion cycle, i.e. 
whenever a total of four successive pulses of the TDC 
signal are generated. This also applies to discrete 
data referred to hereinafter. 

The dead time of the detected air- fuel ratio KACT 
with respect to the target air-fuel ratio KCMD is 
estimated to correspond to about three combustion 
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cycles, and therefore, there is a relationship of 
KCMD(n) = KACT (n+3). When this relationship is 
applied to the equation (23), there is derived an 
equation (24) in FIG. 10. 

Further, the vector 9 (n) of model parameters 
b0(n), rl(n), r2(n), r3(n), and sO(n) in the equation 
(23) is identified with an identification algorithm 
expressed by equations (25) to (31) in FIG. 10. The 
symbol KT(n) in the equation (25) represents a vector 
of a gain coefficient, and the symbol ide_st(n) 
represents an identification error. Further, the 
symbol 0(n) T in the equation (26) represents a 
transposed matrix of 9 (n) . 

The identification error ide_st(n) in the 
equation (25) is calculated using the equations (27) to 
(29) in FIG. 10, and the symbol KACT_HAT(n) in the 
equation (28) represents an identified value of the 
detected air-fuel ratio KACT. Further, the vector KT 
(n) of the gain coefficient is calculated using the 
equation (30) in FIG. 10 /and the symbol T(n) in the 
equation (30) is a square matrix of order 5 defined by 
the equation (31) in FIG. 10. 

In the control system of the present embodiment, 
when the air-fuel ratio control is executed with the 
algorithm expressed by the equations (24) to (31) 
described above, if the LAF sensor 14 has a strong low- 
pass characteristic, the control system can suffer from 
resonance at a repetition period of an integral 
multiple of a period of updating the model parameter 
vector 9 . To overcome the problem, in the present 
embodiment, the second air-fuel ratio controller 40 
calculates the feedback correction coefficient KSTR in 
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the following manner: 

In the second air-fuel ratio controller 40 of the 
control system according to the present embodiment, the 
model parameter vector 0 of the first cylinder #1 
identified by the onboard identifier 41 is oversampled 
in timing synchronous with generation of the TDC signal, 
and at the same time, a moving average value 9 _ave of 
the model parameter vector 9 is calculated. More 
specifically, the moving average value 0_ave(k) of the 
model parameter vector 9 is calculated using an 
equation (32) in FIG. 11, and the feedback correction 
coefficient KSTR(k) is calculated using the moving 
average value 0_ave(k) by an equation (34) in FIG. 11. 
It should be noted that the symbol 6 buf in the 
equation (32) indicates an oversampling value of the 
model parameter vector 9 for the first cylinder #1, 
and the moving average value 6>__ave(k) is defined by an 
equation (33) in FIG. 11. Further, the symbol m in the 
equation (32) represents a predetermined integer, and 
in the present embodiment, m is set to 11. 

As described hereinbefore, the discrete data with 
the symbol (k) in these equations (32) to (34) are data 
sampled in synchronism with the generation of each 
pulse of the TDC signal, and therefore, the 
relationship of n-f = K-4 • f (f: integer) holds. When 
this relationship is applied to the equation (24) in 
FIG. 10, the above equation (34) is derived. Further, 
the identification algorithm with which the model 
parameter vector 9 (k) is identified is expressed by 
equations (35) to (41) shown in FIG. 11. 

As described above, in second air-fuel ratio 
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controller 40 of the control system according to the 
present embodiment, the onboard identifier 40 
identifies the model parameter vector 0 with the 
identification algorithm expressed by the equations 
(35) to (41) in FIG. 11, and the STR controller 42 
calculates the feedback correction coefficient KSTR(k) 
using the equations (32) to (34) in FIG. 11. 

In the following, an air-fuel ratio control 
process, which is executed by the ECU 2, will be 
described with reference to FIGS. 12 to 17. In the 
following description, the symbols (k) and (n) 
indicating that associated values are the current 
values are omitted when deemed appropriate. FIG. 12 
shows a main routine for carrying out the control 
process, which is executed by an interrupt handling 
routine in synchronism with inputting of each pulse of 
the TDC signal. In this process, the final fuel 
injection amount TOUT ± is calculated, on a cylinder-by- 
cylinder basis, as will be described hereinafter. 

First, in a step 1 (in FIG. 12, abbreviated to 
SI; this rule also applies to the other steps referred 
to hereinafter), outputs from the sensors 9 to 19 
described hereinbefore are read in and stored in the 
RAM. 

Then, the process proceeds to a step 2, wherein 
the basic fuel injection amount TIBS is calculated. In 
this process, the basic fuel injection amount TIBS is 
calculated by searching a table, not shown, according 
to the amount of intake air (hereinafter also referred 
to as "the intake air amount") GAIR. 

Then, the process proceeds to a step 3, wherein a 
total correction coefficient KTOTAL is calculated. The 
total correction coefficient KTOTAL is obtained by 
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calculating various correction coefficients by 
searching tables and maps according to various 
operating parameters (e.g. the intake air temperature 
TA, the atmospheric pressure PA, the engine coolant 
temperature TW, the accelerator pedal opening AP, 
throttle valve opening TH, and so forth) and then 
multiplying the thus calculated correction coefficients 
by each other. 

Next, the process proceeds to a step 4, wherein 
the target air-fuel ratio KCMD is calculated. The 
process for calculation of the target air-fuel ratio 
KCMD is not shown here, but is executed by the same 
control method as described in Japanese Laid-Open 
Patent Publication (Kokai) No. 2000-179385. That is, 
the target air-fuel ratio KCMD is calculated depending 
on the operating conditions of the engine 3, by a 
sliding mode control process or a map retrieval process 
such that the output Vout from the 02 sensor 15 
converges to a predetermined target value Vop. 

Then, the process proceeds to a step 5, wherein 
the corrected target air-fuel ratio KCMDM is calculated. 
The corrected target air-fuel ratio KCMDM is calculated 
for compensating for a change in charging efficiency 
due to a change in the air-fuel ratio A/F. The 
corrected target air-fuel ratio KCMDM is calculated by 
searching a table, not shown, according to the target 
air-fuel ratio KCMD calculated in the step 4. 

Next, in steps 6 and 7, the model parameter 
vector 9 of the first cylinder #1 and the feedback 
correction coefficient KSTR are calculated, 
respectively. Processes for calculating these 
parameters will be described in detail hereinafter. 

In the following steps 8 and 10, the vector <f> of 
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the air-fuel ratio variation coefficient, the air-fuel 
ratio variation correction coefficient KOBSVi and 
learned correction value KOBSV_LSi thereof are 
calculated, respectively. Processes for calculating 
these parameters will be described in detail 
hereinafter. 

Then, the process proceeds to a step 11 , wherein 
the demanded fuel injection amount TCYL ± is calculated 
using the basic fuel injection amount TIBS, the total 
correction coefficient KTOTAL, the corrected target 
air-fuel ratio KCMDM, the feedback correction 
coefficient KSTR, the air-fuel ratio variation 
correction coefficient KOBSVi, and the learned 
correction value KOBSV_LS ± thereof, by the following 
equation ( 42 ) : 

TCYL ± = TIBS ■ KTOTAL ' KCMDM* KSTR ' KOBSV ± • KOBSV_LS ± 

.... ( 42 ) 

Then, the process proceeds to a step 12 , wherein 
the final fuel injection amount TOUTi is calculated by 
subjecting the demanded fuel injection amount TCYLi to 
the fuel attachment -dependent correction. More 
specifically, the final fuel injection amount TOUTi is 
calculated by calculating a ratio of an amount of fuel 
attached to the inner walls of the combustion chambers 
to the whole amount of fuel injected from the injectors 
6 during the current combustion cycle, etc. and 
correcting the demanded fuel injection amount TCYLi 
based on the ratio thus calculated. 

Then, the process proceeds to a step 13, wherein 
the drive signal based on the final fuel injection 
amount TOUTi calculated as described above is delivered 
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to the injector 6 for the associated cylinder for which 
the calculation has been performed , followed by 
terminating the present process. 

Next, the process for calculating the model 
parameter vector d executed in the step 6 will be 
described. In this process, first, in a step 20, there 
is carried out a process for setting the cylinder 
number value i which corresponds to the subscript n t " 
in each parameter. 

In this process, the cylinder number value i is 
set based on the immediately preceding value PRVi 
thereof set in the immediately preceding loop as 
follows: When PRVi = 1 holds, the cylinder number 
value i set to 3, when PRVi = 2 holds, the same is set 
to 1, when PRVi = 3 holds, the same is set to 4, and 
when PRVi = 4 holds, the same is set to 2 . As 
described above, the cylinder number value i is 
cyclically set, e.g. in the order of 1 — * 3 — > 4^2~ > 
1^3^4-^2^1 ... 

Next, the process proceeds to a step 21, wherein 
it is determined whether or not the cylinder number 
value i set in the step 20 is equal to a value of 1. 
If the answer to this determination is affirmative 
(YES), which means that the model parameter vector 0 
of the first cylinder #1 is to be calculated, the 
process proceeds to a step 22, wherein the value of the 
model parameter vector 6 calculated and stored in the 
RAM in the immediately preceding loop is set to the 
immediately preceding value PRV 0[6(n-l)]. 

Then, the process proceeds to a step 23, wherein 
the vector £ is calculated using the equation (39) in 
FIG. 11, referred to hereinbefore, and then in a step 
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24, the identified value KACT_HAT of the detected air- 
fuel ratio KACT is calculated using the equation (38) 
in FIG. 11, referred to hereinbefore. 

Then, the process proceeds to a step 25, wherein 
the identification error ide__st is calculated using the 
equation (37) in FIG. 11, referred to hereinbefore, and 
then in a step 26, the value of the next value NEXTT 
[=T(n+l)] of the square matrix calculated and stored 
in the RAM in the immediately preceding loop is set to 
the present value V . 

Then, the process proceeds to a step 27, wherein 
the vector Kf of the gain coefficient is calculated 
using the equation (40) in FIG. 11, referred to 
hereinbefore. Then, the process proceeds to a step 28, 
wherein the model parameter vector 9 is calculated 
using the equation (35) in FIG. 11. 

Then, the process proceeds to a step 29, wherein 
the next value NEXTT of the square matrix is 
calculated by the equation (41) in FIG. 11, and then 
the process proceeds to a step 30, wherein a 
predetermined number (twelve, in the present 
embodiment) of values of the detected air-fuel ratio 
KACT calculated on and before the immediately preceding 
occasion, stored in the RAM, are updated. More 
specifically, each value of the detected air-fuel ratio 
KACT stored in the RAM is set to an older value by one 
control cycle. For example, the current value KACT(k) 
is set to the immediately preceding value KACT(k-l), 
and the immediately preceding value KACT(k-l) is set to 
the second preceding value KACT(k-2), and so forth. 

Then, the process proceeds to a step 31, wherein 
a predetermined number (twelve, in the present 
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embodiment) of oversampling values 9 buf of the model 
parameter vector 9 of the first cylinder #1, stored in 
the RAM, are updated. More specifically, similarly to 
the step 30, each of the oversampling values 9 buf 
stored in the RAM is set to an older value by one 
control cycle. For example, the current oversampling 
value 6>buf(k) is set to the immediately preceding 
oversampling value #buf(k-l), and the immediately 
preceding oversampling value 0buf(k-l) is set to the 
second preceding oversampling value (9buf(k-2), and so 
forth, followed by terminating the present process. 

On the other hand, if the answer to the question 
of the step 21 is negative (NO), which means that it is 
not necessary to calculate the model parameter vector 
9 , the steps 22 to 29 are skipped over, and the steps 
30 and 31 are executed, followed by terminating the 
present process. 

Next, the process for calculating the feedback 
correction coefficient KSTR in the step 7 will be 
described with reference to FIG. 14. In this process, 
first, in a step 40, the moving average value 6>_ave of 
the model parameter vector is calculated based on the 
oversampling values 9 buf updated in the step 31, using 
the equation (32) in FIG. 11. 

Then, in a step 41, the feedback correction 
coefficient KSTR is calculated based on the moving 
average value 0_ave calculated in the step 41, by the 
equation (34) in FIG. 11, referred to hereinbefore. 

Then, the process proceeds to a step 42, wherein 
a predetermined number (twelve in the present 
embodiment) of values of the feedback correction 
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coefficient KSTR calculated in the preceding loops, 
which are stored in the RAM , are updated. More 
specifically, each of the KSTR values stored in the RAM 
is set to an older value by one control cycle. For 
example, the current value KSTR(k) is set to the 
immediately preceding value KSTR(k-l), the immediately 
preceding value KSTR(k-l) is set to the second 
preceding value KSTR(k-2), and so forth. Then, the 
present process is terminated. 

Next, the process for calculating the vector (t> 
of the air-fuel ratio variation coefficient in the step 
8 will be described with reference to FIG. 15. In this 
process, first, in a step 50, the vector <£> of the air- 
fuel ratio variation coefficient calculated in the 
immediately preceding loop, which is stored in the RAM, 
is set to the immediately preceding value PRV<t> [= (f> (k- 
1 ) ] thereof . 

Then the process proceeds to a step 51, wherein 
the vector £ of the simulation values is calculated by 
the equation (7) in FIG. 4, and then to a step 52, 
wherein the estimation value KACT_EST of the detected 
air-fuel ratio is calculated by the equation (6) in FIG. 
4 . 

Then, the process proceeds to a step 53, wherein 
the identification error ide is calculated by the 
equation (5) in FIG. 4, and then to a step 54, wherein 
the value of the next value NEXP [ =P ( k+1 ) ] of the square 
matrix calculated in the immediately preceding loop is 
set to the current value P thereof. 

Then, the process proceeds to a step 55, wherein 
the vector KP of the gain coefficient is calculated by 
the equation (8) in FIG. 4, and then to a step 56, 
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wherein the vector <t> of the air-fuel ratio variation 
coefficient is calculated by the equation (3) in FIG. 4. 

Next, the process proceeds to a step 57, wherein 
the next value NEXP [ =P ( k+1 ) ] of the square matrix is 
calculated by the equation (9) in FIG. 4, and then to a 
step 58, wherein a predetermined number of (12 x 4 in 
the present embodiment) pieces of time-series data of 
the simulation values KACT_OS ± stored in the RAM are 
updated. More specifically, each value of the 
simulation values KACT_OSi stored in the RAM is set to 
an older value by one control cycle (e.g. the current 
value KACT__OSi(k) to the immediately preceding value 
KACT__OSi(k-l ) , and the immediately preceding value 
KACT_OSi(k-l ) to the second preceding value KACT_OS ± (k- 
2), and so forth.) 

Then, the process proceeds to a step 59, wherein 
the current value KACT_OS ± of the simulation value is 
calculated, followed by terminating the present process. 

Next, a process for calculating air-fuel ratio 
variation correction coefficient KOBSV ± executed in the 
step 9 will be described with reference to FIG. 16. In 
this process, first, in the step 70, the moving average 
value Oave of the air-fuel ratio variation coefficient 
is calculated by the equation (10) in FIG. 7. 

Then, the process proceeds to a step 71, wherein 
the following error e is calculated using the equation 
(12) in FIG. 7, referred to hereinbefore, and then in a 
step 72, the integral value Se of the following error 
e is calculated. Then, the process proceeds to a step 
73, wherein the air-fuel ratio variation correction 
coefficient KOBSVi is calculated by the equation (11) 
in FIG. 7, referred to hereinbefore, using the moving 
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average value $ave of the air- fuel ratio variation 
coefficient and the integral value Ze of the following 
error e, calculated in the steps 70 and 72, 
respectively, followed by terminating the present 
process . 

Next, the process for calculating the learned 
correction value KOBSV_LSi of the air-fuel ratio 
variation correction coefficient in the step 10 will be 
described with reference to FIG. 17. In this process, 
first, in a step 80, the exhaust gas volume ESV is 
calculated using the equation (13) in FIG. 9, referred 
to hereinbefore. 

Then, the process proceeds to a step 81, wherein 
the value of the regression coefficient 9 OBSV_LS ± 
calculated in the preceding loop is set to the 
immediately preceding value PRV 0 OBSV_LSi [= 9 
OBSV_LSi(k-l) ] thereof. 

Then, the process proceeds to a step 82, wherein 
the learned correction value KOBSV_LSi is calculated 
using the equation (22) in FIG. 9, referred to 
hereinbefore. Thereafter, the process proceeds to a 
step 83, wherein it is determined whether or not the 
following five conditions (al) to (a5) are satisfied: 

(al) The engine coolant temperature TW is higher 
than a predetermined lower limit value TWOBSL and at 
the same time lower than a predetermined higher limit 
value TWOBSH. 

(a2) The intake air temperature TA is higher 
than a predetermined lower limit value TAOBSL and at 
the same time lower than a predetermined higher limit 
value TWOBSH. 

(a3) The engine speed NE is higher than a 
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predetermined lower limit value NEOBSL and at the same 
time lower than a predetermined higher limit value 
NEOBSH . 

(a4) The intake pipe absolute pressure PBA is 
higher than a predetermined lower limit value PBOBSL 
and at the same time lower than a predetermined higher 
limit value PBOBSH. 

(a5) The vehicle speed VP is higher than a 
predetermined lower limit value VPOBSL and at the same 
time lower than a predetermined higher limit value 
VPOBSH . 

When all of the five conditions (al) to (a5) are 
satisfied, it is judged that the engine is in an 
operating condition in which the regression coefficient 
vector 9 0BSV_LS ± should be calculated by the 
sequential least-squares method, so that the process 
proceeds to a step 84, wherein a vector Z of the 
exhaust gas volume is calculated using the equation 
(19) in FIG. 9, referred to hereinbefore. 

Then, the process proceeds to a step 85, wherein 
the error Eovi is calculated using the equation (17) in 
FIG. 9, referred to hereinbefore, and then to a step 86, 
wherein a next value NEXTQ ± [= Qi(k+1)] of the square 
matrix calculated in the immediately preceding loop, 
stored in the RAM, is set to the current value Q± 
thereof . 

Then, the process proceeds to a step 87, wherein 
the vector KQ ± of the gain coefficient is calculated 
using the equation (20) in FIG. 9, referred to 
hereinbefore, and then to a step 88, wherein the 
regression coefficient vector 9 OBSV_LSi is calculated 
using the equation (15) in FIG. 9, referred to 
hereinbefore. Then, the process proceeds to a step 89, 
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wherein the next value NEXTQi [= Q±(k+1)] of the square 
matrix is calculated using the equation (21) in FIG. 9, 
referred to hereinbefore. 

On the other hand, when the answer to the 
question of the step 83 is negative (NO), i.e. at least 
one of the above five conditions (al) to (a5) is not 
satisfied, the process proceeds to a step 90, wherein 
the immediately preceding value PRV 0 OBSV_LS ± of the 
regression coefficient vector set in the step 81 is set 
to the current value #OBSV_LSi, followed ^by 
terminating the present process. This causes the value 
calculated by the sequential least -squares method in 
the steps 84 to 89 e.g. in the immediately proceeding 
loop to be used as the immediately preceding value PRV 
6 OBSV_LSi of the regression coefficient vector in the 
step 81 in the next loop. 

Next, the operation of the air-fuel ratio control 
executed by the control system 1 will be described with 
reference to FIGS. 18 and 19. FIG. 18 shows an example 
of operation in the case where the air-fuel ratio 
control is carried out by the control system 1 
according to the present embodiment, more specifically, 
the case where during the control provided by the 
second air-fuel ratio controller 40 such that the 
detected air- fuel ratio KACT becomes equal to a value 
of 1 (equivalent ratio corresponding to the 
stoichiometric air-fuel ratio), the first air-fuel 
ratio controller 30 is started from the stopped state, 
in other words, the first air-fuel ratio controller 30 
starts to calculate the air-fuel ratio variation 
coefficient O i# the air-fuel ratio variation 
correction coefficient KOBSVi and the learned 
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correction value KOBSV_LSi thereof. 

FIG. 19 shows, for comparison, a comparative 
example of operation of the air-fuel ratio control in 
the case where the learned correction value KOBSV_LSi 
is calculated with conventional PID control algorithm 
(algorithm expressed by equations (43) and (44) in FIG. 
20) instead of the I-PD control algorithm expressed by 
the equations (11) and (12). In these two figures, the 
values KACTi to KACT 4 represent respective values of 
the air-fuel ratio (values in terms of the equivalent 
ratio) of exhaust gases which are emitted from the 
first to fourth cylinders #1 to #4 and not mixed yet. 
More specifically, the values of KACTi - 4 are calculated 
based on respective outputs from four LAF sensors (not 
shown) for experiment which are additionally disposed 
in the exhaust manifold 7a at respective locations 
immediately downstream of the exhaust ports of the 
cylinders #1 to #4. 

As shown in FIG. 18, in the example of operation 
of the air-fuel ratio control system according to the 
present embodiment, when the first air-fuel ratio 
controller 30 is in stoppage, the values KACTi to KACT 4 
indicative of the respective air-fuel ratio values of 
exhaust gases emitted from the cylinders are made 
unstable, and the detected air-fuel ratio KACT is 
affected thereby and also made somewhat unstable. 
However, when the first air-fuel ratio controller 30 
starts operation (time tl), with the lapse of some time, 
the values KACTi to KACT 4 all converge to a value of 1 
(equivalent ratio corresponding to the stoichiometric 
air-fuel ratio) and accordingly, the detected air-fuel 
ratio KACT also converges to a value of 1. That is, it 
is understood that variation in air-fuel ratio between 
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the cylinders is properly corrected. Further, it is 
understood that the product KOBSVi • KOBSV_LSi ( i = 1 to 
4) of the air-fuel ratio variation correction 
coefficient and the learned correction value thereof is 
also stable. 

In contrast, in the comparative example shown in 
FIG. 19, the setting time from a time point the first 
air-fuel ratio controller 30 starts operation (time t2) 
to a time point all the values of the KACTi_ 4 converge 
to a value of 1 is longer than the example of operation 
of the control system according to the present 
embodiment, and accordingly, the detected air-fuel 
ratio KACT does not smoothly converge to a value of 1 , 
either. In addition, it is understood that the product 
KOBSVi • KOBSV_LSi (i = 1 to 4) of the air-fuel ratio 
variation correction coefficient and the learned 
correction value thereof is not smoothly made stable, 
either. That is, it is understood that compared with 
the case of using the conventional PID control 
algorithm, the use of the I-PD control algorithm as in 
the present embodiment makes it possible to correct 
variation in air- fuel ratio between the cylinders more 
promptly and more appropriately. The reason for this 
is that the learned correction value KOBSV_LSi can be 
more properly calculated with the I-PD control 
algorithm than with the PID control algorithm, without 
causing overshooting in the behavior of the air-fuel 
ratio variation coefficient <&± being caused to 
converge to the moving average value Oave thereof. 

As described above, according to the control 
system of the present embodiment, the first air-fuel 
ratio controller 30 calculates the air-fuel ratio 
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variation coefficient <P±, and calculates the air-fuel 
ratio variation correction coefficient KOBSV ± and the 
learned correction value KOBSV_LSi thereof such that 
the air-fuel ratio variation coefficient 0 ± converges 
to the moving average value $ave. Further, the second 
air-fuel ratio controller 40 calculates the feedback 
correction coefficient KSTR such that the detected air- 
fuel ratio KACT converges to the target air-fuel ratio 
KCMD. Then, the basic fuel injection amount TIBS is 
corrected based on the calculated feedback correction 
coefficient KSTR, air-fuel ratio variation correction 
coefficient KOBSVi, and learned correction value 
KOBSV__LSi thereof, whereby the final fuel injection 
amount TOUTi is calculated on a cylinder-by-cylinder 
basis . 

The adaptive observer 31 of the first air-fuel 
ratio controller 30 estimates the estimation value 
KACT_EST of the detected air- fuel ratio KACT using the 
model [equation (2)] defined by the estimation value 
KACT_EST, the simulation values KACT_OSi, and the air- 
fuel ratio variation coefficients O i# and further 
identifies the air-fuel ratio variation coefficient 0 ± 
as the model parameter by the sequential least -squares 
method such that the estimation value KACT_EST becomes 
equal to the detected air-fuel ratio KACT. This makes 
it possible to remove (filter off) noise-like 
fluctuating components of the exhaust behavior caused 
by a sudden change in the operating condition of the 
engine 3 from the air-fuel ratio variation coefficient 
<E>i, and thereby calculate the air-fuel ratio variation 
coefficient O t as a value substantially indicative of 
variation in air-fuel ratio between the cylinders. 
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Therefore, the basic fuel injection amount TIBS is 
corrected by the variation correction coefficient 
KOBSVi calculated based on the air-fuel ratio variation 
coefficient <J>i, on a cylinder-by-cylinder basis, which 
makes it possible, differently from the conventional 
control system, even when the dynamic characteristics 
of the controlled object are changed due to changes in 
respective contributions of the cylinders to the 
detected air-fuel ratio KACT, which are caused by 
attachment of fuel in the cylinders, variation in the 
response of the LAF sensor 14, and aging of the LAF 
sensor 14, to calculate the final fuel injection amount 
TOUTi on a cylinder-by-cylinder basis, such that 
variation in air-fuel ratio between the cylinders is 
corrected while causing changes in the dynamic 
characteristics of the controlled object to be 
reflected in the model. As a result, according to the 
first embodiment, even in controlling the air-fuel 
ratio of the mixture supplied to the engine 3 having a 
complicated exhaust system layout, it is possible to 
realize a highly robust air-fuel ratio control having a 
large margin of stability, and thereby maintain an 
excellent emission reduction rate of the catalyst. 

Further, the first air-fuel ratio controller 30 
calculates the air-fuel ratio variation correction 
coefficient KOBSV ± with the I-PD control algorithm, 
which makes it possible to calculate the air- fuel ratio 
variation correction coefficient KOBSVi such that 
overshooting is not caused in the behavior of the air- 
fuel ratio variation coefficient 0 ± being caused to 
converge to the moving average value Oave thereof. 
This makes it possible to correct variation in air-fuel 
ratio between the cylinders, while preventing the air- 
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fuel ratio of exhaust gases from each cylinder from 
exhibiting an oscillatory behavior. Further, since the 
air-fuel ratio variation correction coefficient KOBSVi 
is calculated such that the air-fuel ratio variation 
coefficient 4> ± is caused to converge to the moving 
average value Oave thereof, it is possible to correct 
variation in air- fuel ratio between the cylinders while 
preventing the air-fuel ratio by the first air-fuel 
ratio controller 30 and the air-fuel ratio control by 
the second air-fuel ratio controller 40 from 
interfering with each other. 

Furthermore, the first air-fuel ratio controller 
30 calculates the learned correction value KOBSV_LS ± of 
the air-fuel ratio variation correction coefficient 
KOBSVi by the regression equation [equation (22)] using 
the exhaust gas volume ESV as an independent variable, 
and at the same time, the regression coefficient vector 
<9 0BSV_LSi as a vector of the regression coefficient 
AOBSV_LSi and the constant term BOBSV_LSi is calculated 
by the sequential least-squares method. This makes it 
possible to calculate the learned correction value 
KOBSV_LSi as a value in which the variation in air- fuel 
ratio between the cylinders is properly reflected 
therein, even when the engine 3 is in an operating 
condition which can change drastically, such as a 
transient operating condition, causing a sudden change 
in the state of variation in air-fuel ratio between the 
cylinders. Therefore, even when the engine 3 is in a 
transient operating condition, it is possible to 
control the air-fuel ratio while compensating for the 
variation in air-fuel ratio between the cylinders. 

Further, the air-fuel ratio variation coefficient 
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<£>i and the regression coefficient vector #OBSV__LSi are 
calculated by the sequential least-squares method, and 
therefore, compared with the case of using the general 
least-squares method as the statistical algorithm, it 
is possible to calculate the air-fuel ratio variation 
correction coefficient KOBSV ± and the learned 
correction value KOBSV_LSi every control cycle, even at 
the start of the air-fuel ratio control. Therefore, by 
setting the initial values of the air-fuel ratio 
variation correction coefficient KOBSVi and the learned 
correction value KOBSV_LS± in advance, it is possible 
to calculate the final fuel injection amount TOUT ± as a 
value always corrected by the product of the learned 
correction value KOBSV_LS ± calculated every control 
cycle and the air-fuel ratio variation correction 
coefficient KOBSVi at the start of the air-fuel ratio 
control, whereby the controllability at the start of 
the air-fuel ratio control can be enhanced. This makes 
it possible to enhance the emission reduction rate of 
the catalyst at the start of the air-fuel ratio control. 

Although in the first embodiment, the first- 
degree equation is employed as the regression equation 
used in the calculation of the learned correction value 
KOBSV__LSi, this is not limitative, but an n-th-degree 
equation (n is an integer equal to or larger than 2) 
may be used. In such a case as well, by calculating 
the regression coefficients and the constant terms of 
the n-th-degree equation by the sequential least - 
squares method, it is possible to obtain the same 
advantageous effects as provided by the first 
embodiment. Further, the learned correction value 
KOBSV_LSi may be calculated by using predetermined 
values set for each of a plurality of operating regions 
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in advance, as the regression coefficient and the 
constant term of the regression equation. This can 
reduce the time for computing the learned correction 
value KOBSV_LSi, thereby reducing the computational 
load on the ECU 2 . 

Further, although in the first embodiment, the I- 
PD control algorithm is employed as the control 
algorithm for causing the air-fuel ratio variation 
coefficient $i to converge to the moving average value 
Oave thereof, this is not limitative, but it goes 
without saying that another suitable algorithm may be 
employed. For example, instead of the I-PD control 
algorithm, an IP-D control algorithm (differential- 
preceding PID control algorithm) expressed by equations 
(45) and (46) in FIG. 20 may be employed to calculate 
the air-fuel ratio variation correction coefficient 
KOBSVi or a response-specified control algorithm 
(sliding mode control algorithm or back- stepping 
control algorithm) expressed by equations (47) to (49) 
in FIG. 20 may be employed to calculate the air- fuel 
ratio variation correction coefficient KOBSVi. Even 
when one of these control algorithms is employed, 
similarly to the case of the present embodiment using 
the I-PD control algorithm, it is possible to calculate 
the air-fuel ratio variation correction coefficient 
KOBSVi such that overshooting is not caused in the 
behavior of the air-fuel ratio variation coefficient O 
i being caused to converge to the moving average value 
$ave thereof. As a result, it is possible to promptly 
and appropriately correct the variation in air-fuel 
ratio between the cylinders. 

Further, as described above, when the I-PD 



99 



H02-2749US 



control algorithm, IP-D control algorithm, and the 
response-specified control algorithm are employed in 
the calculation of the air-fuel ratio variation 
correction coefficient KOBSV i# the feedback gain 
thereof may be determined based on the optimal 
regulator theory or the Hoo control theory. This makes 
it possible to more effectively suppress overshooting 
in the behavior of the air-fuel ratio variation 
coefficient 0 ± being caused to converge to the moving 
average value Oave thereof, with the result that the 
accuracy of correction of variation in air-fuel ratio 
between the cylinders can be further enhanced. 

Further, it goes without saying that when the 
setting time over which the air-fuel ratio variation 

coefficient 0± converges to the moving average value <J> 
ave thereof may be long, the air- fuel ratio variation 
correction coefficient KOBSVi may be calculated with 
the PID control algorithm expressed by the equations 
(43) and (44) in FIG. 20. Further, the average value 
of the variation coefficient as a target value to which 
the air-fuel ratio variation coefficient $i is caused 
to converge is not limited to the moving average value 
$ave in the present embodiment, but it may be a 
weighted average value. 

Further, in the illustrated example of the first 
embodiment, the adaptive observer 31 of the first air- 
fuel ratio controller 30 identifies the vector <Mk) of 
the air-fuel ratio variation coefficient, by the 
variable-gain sequential least-squares method expressed 
by the equations (3) to (9) in FIG. 4 shown in FIG. 4, 
it goes without saying that the identification 
algorithm with which the adaptive observer 31 
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identifies the vector (f> (k) of the air-fuel ratio 
variation coefficient is not limited to this. For 
example, the air-fuel ratio variation coefficient may 
be identified by the fixed gain method to which is 
applied the 6 correction method expressed by equations 
(50) to (57) in FIG. 21. 

The symbol $base in the equation (50) in FIG. 21 
represents a reference value vector (model parameter 
reference value) defined by the equation (51), and four 
elements of this vector, i.e. reference values ( Pbasel 
to <I ) base4 are calculated by searching a table shown in 
FIG. 2 2 according to the exhaust gas volume ESV. As 
shown in FIG. 22, the four reference values <E>basel to 
Obase4 are all set to a value close to a value of 1. 
Further, the symbol d(J>(k) in the equation (50) 
represents a correction term (correction component) 
defined by the equation (52), and calculated by the 
equations (53) to (57). 

When the vector (£> (k) of the air-fuel ratio 
variation coefficient is identified by the fixed gain 
method to which the 6 correction method is applied, 
the computing time can be reduced compared with the 
case of using the sequential least-squares method, and 
the computational load on the ECU 2 can be reduced. As 
a result, it is possible to reduce the size and cost of 
the ECU 2. Moreover, even when the engine 3 is in an 
operating condition in which the air-fuel ratio changes 
violently, such as a transient operating condition, the 
identified value of the vector <Mk) can be constrained 
to values close to a value of 1 , which makes it 
possible to promptly and properly calculate the vector 
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(/>(k) of the air-fuel ratio variation coefficient 
representative of the variation in air-fuel ratio 
between the cylinders as a value in which the behavior 
of the air-fuel ratio is properly reflected, whereby 
the stability of the air-fuel ratio control can be 
enhanced. 

When the tables shown in FIG. 2 2 cannot be 
provided in advance, all the four elements Obasel to O 
base4 may be set to a value of 1. 

Further, although in the present embodiment, the 
basic fuel injection amount TIBS is calculated by 
searching the table according to the intake air amount 
GAIR, in the step 2 in FIG. 12, this is not limitative, 
but the basic fuel injection amount TIBS may be 
calculated by searching a map according to the intake 
pipe absolute pressure PBS and the engine speed NE . 

Next, a description will be given of the control 
system 101 according to a second embodiment of the 
present invention. The control system 101 is 
distinguished from the control system 101 according to 
the first embodiment in that as shown in FIGS. 23 and 
24 by a third air-fuel ratio controller 60 which is 
additionally provided, and the remainder is identical 
in construction. Therefore, the following description 
will be mainly given of the third air-fuel ratio 
controller 60 (third fuel amount -correcting means), 
with the component parts identical to those of the 
first embodiment being designated by the same reference 
numerals, and description thereof being omitted unless 
otherwise required . 

In the control system 101, as described 
hereinafter, the third air-fuel ratio controller 60 
calculates an intake air amount variation correction 
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coefficient KICYLi and a learned correction value 
KICYL_LS ± thereof so as to correct variation in intake 
air amount between the cylinders. Then, the basic fuel 
injection amount TIBS is multiplied by the corrected 
target air-fuel ratio KCMDM , the total correction 
coefficient KTOTAL , the feedback correction coefficient 
KSTR , the air-fuel ratio variation correction 
coefficient KOBSVi, the learned correction value 
KOBSV_LSi of the air-fuel ratio variation correction 
coefficient, the intake air amount variation correction 
coefficient KICYL i# and a learned correction value 
KICYL_LS± of the intake air amount variation correction 
coefficient, whereby the demanded fuel injection amount 
TCYLi is calculated, on a cylinder-by-cylinder basis. 
Then, the fuel attachment-dependent correction section 
50 calculates the final fuel injection amount TOUTi 
based on the demanded fuel injection amount TCYLi, on a 
cylinder-by-cylinder basis. 

Next, a description will be given of the third 
air-fuel ratio controller 60. As shown in FIG. 25, 
when the intake air amount GAIR to be supplied to the 
engine 3 is detected by the air flow sensor 9, 
pulsation of the intake air caused by the suction 
behavior of each cylinder is also detected. When there 
occurs variation in intake air amount between the 
cylinders, the pulsation of intake air becomes 
irregular as shown in FIG. 25. That is, FIG. 25 shows 
a case in which the intake air amount in the fourth 
cylinder #4 is smaller than those of the other 
cylinders . 

This air-fuel ratio controller 60 estimates the 
variation in intake air between the cylinders , for 
correction of the fuel injection amount based on the 
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estimated variation, and is comprised of an adaptive 
observer 61, an intake air amount variation correction 
coefficient-calculating section 62, a learned 
correction value-calculating section 63, and a 
multiplication section 64. In this third air-fuel 
ratio controller 60, with algorithms described 
hereinbelow, the adaptive observer 61 (simulation 
value -generating means, estimation means, 
identification means, delay means) calculates an intake 
air amount variation coefficient M^i, on a cylinder-by- 
cylinder basis, and the intake air amount variation 
correction coefficient-calculating section 62 (first 
control means, third correction value-calculating 
means) calculates the intake air amount variation 
correction coefficient KICYL i# on a cylinder-by- 
cylinder basis. Further, the learned correction value- 
calculating section 63 (learned correction value- 
calculating means) calculates the learned correction 
value KICYL_LSi of the intake air amount variation 
correction coefficient, on a cylinder-by-cylinder basis. 
Further, the multiplication section 64 (correction 
means) multiplies the intake air amount variation 
correction coefficients KICYLi to KICYL 4 by the learned 
correction values KICYL_LSi to KICYL_LS 4 , respectively, 
that is, the intake air amount variation correction 
coefficient KICYLi is corrected by the learned 
correction value KICYL_LS ± . 

Next, a description will be given of the 
algorithm of the adaptive observer 61. First, as shown 
in FIG. 26, the intake system of the engine 3 is 
regarded as a system which is represented by four 
simulation values GAIR_0Si to GAIR_OS 4 and four intake 
amount variation coefficients to These 
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simulation values GAIR_OSi to GAIR_OS 4 are values 
simulating the intake start timing of intake air and 
the intake air behavior, on a cylinder-by-cylinder 
basis , and the intake air amount variation coefficient 

represents variation in intake air amount between 
the cylinders and the amount of change in the intake 
air behavior. When this system is modeled into a 
discrete- time system model, there is obtained an 
equation (58) in FIG. 27. In the equation (58), the 
symbol d' represents dead time (predetermined delay 
time) which the air takes to flow in the intake pipe 4 
from the air flow sensor 9 to each cylinder, and is set 
to a predetermined fixed value in the present 
embodiment. The dead time d' may be set depending on 
an operating condition (e.g. the engine speed NE ) of 
the engine 3 . 

The adaptive observer 61 according to the present 
embodiment uses an equation formed by replacing the 
left side of the equation (58) by the estimation value 
GAIR_EST(k) of the intake air amount, i.e. a model 
represented by an equation (59) in FIG. 27, and a 
signal generator 61a generates the simulation value 
GAIR_OSi, as described hereinafter. At the same time, 
a vector ip (k) of the intake air amount variation 
coefficient W ± as a model parameter of the equation 
(59) is identified by the variable-gain sequential 
least -squares method expressed by equations (60) to 
(66) in FIG. 27 such that the estimation value 
GAIR_EST(k) becomes equal to the intake air amount 
GAIR(k-d' ) . 

The symbol KP(k) in the equation (60) represents 
a vector of a gain coefficient, and the symbol ide'(k) 
represents an identification error. Further, the 
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symbol if) (k) T in the equation (61) represents a 
transposed matrix of tl) (k) . The identification error 
ide'(k) in the equation (60) is calculated by equations 
(62) to (64) in FIG. 27, and the symbol £' (k) in the 
equation (63) represents a vector of the simulation 
value defined by an equation (64). Further, the vector 
KP(k) of the gain coefficient is calculated by an 
equation (65) in FIG. 27, and the symbol R(k) in the 
same equation is a square matrix of order 4 defined by 
an equation (66) in FIG. 27. 

As described above, this adaptive observer 61 
identifies the vector ii> (k) of the intake air amount 
variation coefficient ?i with the algorithm based on 
the sequential least -squares method shown in the 
expressions (60) to (66). This makes it possible to 
remove (filter off) noise-like fluctuating components 
of the intake air behavior caused by a sudden change in 
the operating condition of the engine 3 from the intake 
air amount variation coefficient ? 1# and thereby 
calculate the intake air amount variation coefficient 
M^i as a value substantially indicative of variation in 
intake air amount between the cylinders . 

The configuration of the adaptive observer 61 can 
be represented by a block diagram shown in FIG. 28, 
similarly to the adaptive observer 31 of the first air- 
fuel ratio controller 30. That is, as shown in FIG. 28, 
in the adaptive observer 61, the signal generator 61a 
generates the vector £ ' (k) of the simulation values 
GAIR_OSi. More specifically, as shown in FIG. 29, the 
signal generator 61a generates the simulation value 
GAIR_OSi as a signal values such having a waveform of a 
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combination of triangular waves and trapezoidal waves 
formed such that the total sum of the simulation values 
constantly becomes equal to a value of 1. Further, the 
multiplier 61b generates the estimation value 
GAIR_EST(k) of the intake air amount as a value 
obtained by multiplying the vector £ ' (k) of the 
simulation values by the vector ^(k-1) of the intake 
air amount variation coefficient. Then, the 
differentiator 61d generate the identification error 
ide'(k) as the difference between the intake air amount 
GAIR(k-d') and the estimation value GAIR__EST(k ) . 

Further, a logic unit 61e generates the vector 
KP(k) of the gain coefficient based on the vector 
£ '(k) of the simulation values, and a multiplier 61f 
generates the product [ide'(k) • KP(k) ] of the 
identification error ide'(k) and the vector KP(k) of 
the gain coefficient. Next, an adder 61g generates the 
vector il> (k) of the intake air amount variation 
coefficient as the sum of the product [ide'(k) 'KP(k)] 
and the delayed vector ^(k-1) of the intake air amount 
variation coefficient . 

Next, an algorithm with which the intake air 
amount variation correction coefficient -calculating 
section 62 calculates the intake air amount variation 
correction coefficient KICYLi( first input, third 
correction value). In the intake air amount variation 
correction coefficient-calculating section 62, first, 
by an equation (67) in FIG. 30, the moving average 
value ?ave(k) of the intake air amount variation 
coefficient is calculated based on the vector ?/>(k) of 
the intake air amount variation coefficient calculated 
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by adaptive observer 61, i.e. the four intake air 
amount variation coefficients ^(k) to ^ 4 (k). Next, 
the intake air amount variation correction coefficient 
KICYLi is calculated by the I-PD control 
(proportional/differential-preceding PID control) 
algorithm, on a cylinder-by-cylinder basis, such that 
the intake air amount variation coefficient ^i(k) is 

caused to converge to the moving average value ?ave(k) . 
This I-PD control algorithm is expressed by 
equations ( 68 ) and (69) in FIG. 30. 

Here, the air-fuel ratio control for correcting 
variation in intake air amount between the cylinders by 
the third air-fuel ratio controller 60 has a 
possibility of interfering with the air-fuel ratio 
control for correcting the variation in air-fuel ratio 
between the cylinders by the first air-fuel ratio 
controller 30. To avoid this inconvenience, it is 
necessary to make the speed at which the controller 30 

causes the intake air amount variation coefficient W 
i(k) to converge to the moving average value Oave and 
the speed at which the controller 60 causes the intake 
air amount variation coefficient W± (k) to converge to 

the moving average value ?ave(k) , different from each 
other . 

In the present embodiment, the feedback gains FI ' , 
GI ' , and HI' in the above equation (68) are set such 
that the absolute values thereof are larger than the 
corresponding absolute values of the feedback gains FI , 
GI, HI in the equation (11) referred to hereinabove. 
In other words, the feedback gains FI ' , GI ' , and HI' 
are set such that there is satisfied the relationship 
of 0 < |FI| < |FI'|, 0 < | GI | < |GI'|, and 0 < | HI | < 
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|HI'|. This makes it possible to control the air-fuel 
ratio such that the converging speed at which the 
intake air amount variation coefficient ^i(k) is 
caused to converge to the moving average value ?ave(k) 
is faster than the converging speed at which the air 
fuel ratio variation coefficient Oi(k) is caused to 
converge to the moving average value Oave. This is 
because due to a higher S/N ratio of the air flow 
sensor 9 than that of the LAF sensor 14, by setting the 
feedback gains to satisfy the above-mentioned 
relationship, the stability of the air-fuel ratio 
control can be ensured as a whole while preventing the 
two types of the air- fuel ratio control from 
interfering with each other. 

Moreover, the intake air amount variation 
correction coefficient KICYLi is calculated with the I- 
PD control algorithm such that the intake air amount 
variation coefficient ^i(k) is caused to converge to 
the moving average value ?ave(k) thereof. This makes 
it possible to provide control such that overshooting 
is not caused in the behavior of the intake air amount 
variation coefficient ^±(k) being caused to converge 
to the moving average value ?ave(k) . This makes it 
possible to prevent the controllability of the second 
air-fuel ratio controller 40 for the air-fuel ratio 
control from being lowered, when the third air- fuel 
ratio controller 60 carries out the air-fuel ratio 
control for correcting variation in intake air amount 
between the cylinders . 

Next, a description will be given of an algorithm 
with which the learned correction value -calculating 
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section 63 calculates the learned correction value 
KOBSV_LSi of the intake air amount variation correction 
coefficient KICYL ± . The intake air amount variation 
correction coefficient KICYLi is susceptible to 
operating conditions of the engine 3, and when the 
operating conditions of the engine 3 are changed, the 
coefficient KICYLi is changed accordingly. FIG. 31 
shows the relationship between the exhaust gas volume 
ESV(k) as an operating condition parameter indicative 
of an operating condition of the engine and the intake 
air amount variation correction coefficient KICYLi(k) . 

Referring to FIG. 31, it can be seen, similarly 
to the air-fuel ratio variation correction coefficient 
KOBSVi(k), that an approximate value i.e. estimation 
value of the intake air amount variation correction 
coefficient KICYLi(k) can be calculated by a first- 
degree equation using the intake air amount variation 
correction coefficient KICYL ± (k) as a dependent 
variable and the exhaust gas volume ESV(k) as an 
independent variable. Therefore, in the learned 
correction value-calculating section 63, the learned 
correction value KICYL_LSi(k) of the intake air amount 
variation correction coefficient is defined as an 
estimation value calculated by a regression equation 
expressed by an equation (70) in FIG. 32, and a vector 
# ICYL_LSi(k) of a regression coefficient AICYL_LSi and 
a constant term BICYL_LS ± (hereinafter referred to as 
"the regression coefficient vector") is calculated by a 
sequential least -squares method expressed by equations 
(71) to (77) shown in FIG. 32. 

In this equation (71), the symbol KUi(k) 
represents a vector of a gain coefficient, and the 
symbol Eic ± (k) represents an error. Further, the error 
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Eici(k) is calculated using an equation (73) in FIG. 32. 
Further, the vector KUi(k) of the gain coefficient is 
calculated using ah equation (76) in FIG. 32, and the 
symbol U ± (k) in this equation (76) represents a square 
matrix of order 2 defined by an equation (77) in FIG. 
32 . 

Further, the learned correction value 
KICYL_LSi(k) is more specifically calculated using an 
equation (78) in FIG. 32. It should be noted that when 
the engine 3 is in an extreme operating condition or 
operating environment, the calculation of the 
regression coefficient AICYL_LSi and the constant term 
BICYL_LSi by the sequential least -squares method is 
avoided, and the immediately preceding value 9 
ICYL_LSi(k- 1 ) of the regression coefficient vector is 
used as the current value 0ICYL_LS ± (k) in calculation 
of the learned correction value KICYL_LS ± (k) . 

With the algorithm expressed by the equations 
(71) to (78), the learned correction value-calculating 
section 63 calculates the learned correction value 
KICYL_LSi(k) such that the learned correction value 
KICYL_LSi(k) converges to the product of the learned 
correction value KICYL_LSi(k) and the intake air amount 
variation correction coefficient KICYLi(k) . 

It should be noted that as shown in FIG. 25, even 
when the intake air absolute pressure PBA is detected 
by the intake pipe absolute pressure sensor 11, it is 
possible to detect pulsation of intake air, so that 
with an algorithm in which the parameter represented by 
the intake air amount GAIR in the above equations (58) 
to (78) is replaced by a parameter represented by the 
intake pipe absolute pressure PBA, and using the intake 
pipe absolute pressure PBA detected by the intake pipe 
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absolute pressure sensor 11, it is possible to form an 
air-fuel ratio controller for correcting variation in 
the intake air amount between the cylinders. 

In the following, an air-fuel ratio control 
process according to the second embodiment will be 
described with reference to FIGS. 33 to 36. FIG. 33 
shows a main routine for carrying out the control 
process, which is executed by an interrupt handling 
routine in synchronism with inputting of each pulse of 
the TDC signal. As shown in FIG. 33, the steps other 
than steps 111 to 113 are identical to the steps 1 to 
13 in FIG. 12, and therefore, a description will be 
given mainly of the steps 111 to 113. 

More specifically, in a step 110, the learned 
correction value KOBSV_LS ± of the air- fuel ratio 
variation correction coefficient is calculated, and 
then, the vector ii> of the intake air amount variation 
coefficient is calculated as described in detail 
hereinafter. 

Then, the process proceeds to a step 112, wherein 
the intake air amount variation correction coefficient 
KICYLi is calculated, and to a step 113, wherein the 
learned correction value KICYL_LSi of the intake air 
amount variation correction coefficient is calculated. 
Next, similarly to the steps 11 to 13, the steps 114 to 
116 are executed, followed by terminating the present 
process . 

Next, the process for calculating the vector ip 
of the intake air amount variation coefficient executed 
in the step 111 will be described in detail with 
reference to FIG. 34. In this process, by the same 
method of calculation of the vector <t> of the air-fuel 
ratio variation coefficient described hereinabove with 
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reference to FIG. 15, the vector ifr of the intake air 
amount variation coefficient is calculated. More 
specifically, in a step 120, the vector ii> of the 
intake air amount variation coefficient calculated in 
the immediately preceding loop and stored in the RAM is 
set to the immediately preceding value PRV^ [= *Mk-l)]. 

Then, the process proceeds to a step 121, wherein 
the current value GAIR_0S ± of the simulation value is 
calculated, and then to a step 122, wherein the vector 
£ of the simulation value is calculated by the 
equation (64) in FIG. 27. 

Then, the process proceeds to a step 123, wherein 
the estimation value GAIR_EST of the intake air amount 
is calculated by the equation (63) in FIG. 27, and then 
to a step 124, wherein the identification error ide ' is 
calculated by the equation (62) in FIG. 27. 

Then, the process proceeds to a step 125, wherein 
the next value NEXR [=R(k+l)] of the square matrix 
calculated in the immediately preceding loop is set to 
the current value R thereof, and then to a step 126, 
wherein the vector KR of the gain coefficient is 
calculated by the equation (65) in FIG. 27. 

Then, the process proceed to a step 127, wherein 
the vector tp of the intake air amount variation 
coefficient is calculated by the equation (60) in FIG. 
27, and then to a step 128, wherein the next value NEXR 
[=R(k+l)] of the square matrix is calculated by the 
equation (66) in FIG. 27. 

Then, the process proceeds to a step 129, wherein 
a predetermined number (twelve, in the present 
embodiment) of values of the intake air amount GAIR 
stored in the RAM, are updated. More specifically, 
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each value of the intake air amount GAIR stored in the 
RAM is set to an older value by one control cycle. For 
example, the current value GAIR(k) is set to the 
immediately preceding value GAIR(k-l), the immediately 
preceding value GAIR(k-l) is set to the second 
preceding value GAIR(k-2), and so forth, followed by 
terminating the present process. 

Next, the process for calculating the intake air 
amount variation correction coefficient KICYL ± in the 
step 112 will be described with reference to FIG. 35. 
In this process, the intake air amount variation 
correction coefficient KICYL ± is calculated in a manner 
similar to the calculation of the air-fuel ratio 
variation correction coefficient KOBSV ± , described 
hereinbefore with reference to FIG. 16. More 
specifically, first, in a step 140, the moving average 
value ?ave of the intake air amount variation 
coefficient is calculated by the equation (67) in FIG. 
30. 

Then, the process proceeds to a step 141, wherein 
the following error e' is calculated by the equation 
(69) in FIG. 30, and then to a step 142, wherein the 
integral value Se' of the following error is 
calculated. Then, the process proceeds to a step 143, 
wherein the air-fuel ratio variation correction 
coefficient KICYL ± is calculated using the moving 
average value fave of the intake air amount variation 
coefficient and the integral value Ze' of the 
following error calculated in the steps 140 and 142, 
respectively, by the equation (68) in FIG. 30, followed 
by terminating the present process. 

Next, the process for calculating the learned 
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correction value KICYL_LSi of the intake air amount 
variation correction coefficient in the step 113 will 
be described with reference to FIG. 36. In this 
process, the learned correction value KICYL_LSi of the 
intake air amount variation correction coefficient is 
calculated in a manner similar to the calculation of 
the learned correction value KOBSV_LSi of the air- fuel 
ratio variation correction coefficient described 
hereinbefore with reference to FIG. 17. More 
specifically, first, in a step 150, the exhaust gas 
volume ESV is calculated using the equation (13) in FIG. 
9, referred to hereinbefore. 

Then, the process proceeds to a step 151, wherein 
the value of the regression coefficient vector 9 
ICYL_LSi calculated in the immediately preceding loop 
is set to the immediately preceding value PRV0ICYL_LSi 
[= 9 ICYL_LS ± (k-l ) ] thereof. 

Then, the process proceeds to a step 152, wherein 
the learned correction value KICYL_LSi is calculated 
using the equation (78) in FIG. 32, referred to 
hereinbefore. Thereafter, the process proceeds to a 
step 153, wherein it is determined whether or not all 
the following five conditions ( a6 ) to (alO) are 
satisfied: 

(a6) The engine coolant temperature TW is higher 
than a predetermined lower limit value TWICYL and at 
the same time lower than a predetermined higher limit 
value TWICYH. 

(a7) The intake air temperature TA is higher 
than a predetermined lower limit value TAICYL and at 
the same time lower than a predetermined higher limit 
value TWICYH. 

(a8) The engine speed NE is higher than a 
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predetermined lower limit value NEICYL and at the same 
time lower than a predetermined higher limit value 
NEICYH . 

(a9) The intake pipe absolute pressure PBA is 
higher than a predetermined lower limit value PBICYL 
and at the same time lower than a predetermined higher 
limit value PBICYH. 

(alO) The vehicle speed VP is higher than a 
predetermined lower limit value VPICYL and at the same 
time lower than a predetermined higher limit value 
VPICYH. 

When all of the five conditions (a6) to (alO) are 
satisfied, it is judged that the engine is in an 
operating condition in which the regression coefficient 
vector 9 ICYL_LSi should be calculated by the 
sequential least -squares method, so that the process 
proceeds to a step 154, wherein a vector Z' of the 
exhaust gas volume is calculated using the equation 
(75) in FIG. 32, referred to hereinbefore. 

Then, the process proceeds to a step 155, wherein 
the error Eic ± is calculated using the equation (73) in 
FIG. 32, referred to hereinbefore, and then to a step 
156, wherein a next value NEXUi [= U ± (k+1)] of the 
square matrix calculated in the immediately preceding 
loop, stored in the RAM, is set to the current value U ± 
thereof . 

Then, the process proceeds to a step 157, wherein 
the vector KU ± of the gain coefficient is calculated 
using the equation (76) in FIG. 32, referred to 
hereinbefore, and then to a step 158, wherein the 
regression coefficient vector 9 ICYL_JLSi is calculated 
using the equation (71) in FIG. 32, referred to 
hereinbefore. Then, the process proceeds to a step 159, 
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wherein the next value NEXU ± [= Ui(k+1)] of the square 
matrix is calculated using the equation (77) in FIG. 32, 
referred to hereinbefore, followed by terminating the 
present process. 

On the other hand, when the answer to the 
question of the step 153 is negative (NO), i.e. at 
least one of the above five conditions (a6) to (alO) is 
not satisfied, the process proceeds to a step 160, 
wherein the immediately preceding value PRV 0 ICYL_LSi 
of the regression coefficient vector set in the step 
151 is set to the current value 0ICYL_LS ± , followed by 
terminating the present process. This causes the value 
calculated by the sequential least - squares method in 
the steps 154 to 159 e.g. in the immediately preceding 
loop to be used as the immediately preceding value PRV 
0 ICYL_LSi of the regression coefficient vector in the 
step 151 in the next loop. 

As described above, according to the control 
system 101 of the second embodiment, the first air-fuel 
ratio controller 30 calculates the air-fuel ratio 
correction coefficient KOBSV ± and the learned 
correction value KOBSV_LSi thereof, and the second air- 
fuel ratio controller 40 calculates the feedback 
correction coefficient KSTR. Further, the third air- 
fuel ratio controller 60 calculates the intake air 
amount variation coefficient V ± and calculates the 
intake air amount variation correction coefficient 
KICYLi and the learned correction value KICYL_LSi 
thereof such that the intake air amount variation 
coefficient ? A converges to the moving average value W 
ave thereof. Then, the basic fuel injection amount 
TIBS is corrected by the calculated feedback correction 



117 



H02-2749US 



coefficient KSTR , the air-fuel ratio variation 
correction coefficient KOBSV i# the learned correction 
value KOBSV_LSi, the intake air amount variation 
correction coefficient KICYL i# and the learned 
correction value KICYL_LS i# whereby the final fuel 
injection amount TOUTi is calculated, on a cylinder-by- 
cylinder basis. 

The adaptive observer 61 of the third air- fuel 
ratio controller 60 estimates the estimation value 
GAIR_EST of the intake air amount GAIR using a model 
[equation (59)] defined by the estimation value 
GAIR__EST , the simulation values GAIR_0S i# and the 
intake air amount variation coefficients ? i# and 
further the intake air amount variation coefficient ?i 
as the model parameter by the sequential least -squares 
method such that the estimation value GAIR_EST becomes 
equal to the intake air amount GAIR. This makes it 
possible to remove (filter off) noise-like fluctuating 
components of the exhaust behavior caused by a sudden 
change in the operating condition of the engine 3 from 
the intake air amount variation coefficient and 
thereby calculate the intake air amount variation 
coefficient W± as a value substantially indicative of 
variation in intake air amount between the cylinders . 
Therefore, since the basic fuel injection amount TIBS 
is corrected by the intake air amount variation 
correction coefficient KICYVi calculated based on the 
intake air amount variation coefficient ^i, on a 
cylinder -by- cylinder basis, it is possible, even when 
the dynamic characteristics of the controlled object 
are changed due to variation in the response of the air 
flow sensor 9 and the aging of the same, to calculate 
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the final fuel injection amount TOUT i# on a cylinder- 
by-cylinder basis, such that variation in intake air 
amount between the cylinders is corrected while causing 
changes in the dynamic characteristics of the 
controlled object to be reflected in the model. As a 
result, according to the present embodiment, even in 
controlling the air-fuel ratio of the mixture supplied 
to the engine 3 having a complicated exhaust system 
layout, it is possible to realize a highly robust air- 
fuel ratio control having a large margin of stability, 
and thereby maintain an excellent emission reduction 
rate of the catalyst. 

Further, the third air-fuel ratio controller 60 
calculates the intake air amount variation correction 
coefficient KICYLi with the I-PD control algorithm, 
which makes it possible to calculate the intake air 
amount variation correction coefficient KICYL ± such 
that overshooting is not caused in the behavior of the 
intake air amount variation coefficient ? A being 
caused to converge to the moving average value Vave 
thereof. This makes it possible to correct variation 
in intake air amount between the cylinders, while 
preventing the amount of intake air supplied to each 
cylinder from exhibiting an oscillatory behavior. 
Moreover, in the I-PD control algorithm, the feedback 
gains FI ' , GI ' , and HI' are set to respective values 
such that the speed at which the intake air amount 
variation coefficient W ± converges to the moving 
average value ^Fave becomes higher than the speed at 
which the air-fuel ratio variation coefficient 0 ± 
converges to the moving average value <E>ave, so that it 
is possible to prevent the air-fuel ratio control by 
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the third air-fuel ratio controller 60 and the air-fuel 
ratio control by the first air-fuel ratio controller 30 
from interfering with each other. Further, since the 
intake air amount variation correction coefficient 
KICYLi is calculated such that the intake air amount 
variation coefficient ?i is caused to converge to the 
moving average value ^ave thereof, it is possible to 
prevent the air-fuel ratio control by the third air- 
fuel ratio controller 60 and the air-fuel ratio control 
by the second air-fuel ratio controller 40 from 
interfering with each other. Thus, the variation in 
intake air amount between the cylinders can be 
corrected while preventing the air-fuel ratio control 
by the third air-fuel ratio controller 60, the air-fuel 
ratio control by the first air-fuel ratio controller 30, 
and the air-fuel ratio control by the second air-fuel 
ratio controller 40 from interfering with each other. 

Furthermore, the third air-fuel ratio controller 
60 calculates the learned correction value KICYL_LS ± of 
the intake air amount variation correction coefficient 
KICYLi by the regression equation [equation (78)] using 
the exhaust gas volume ESV as an independent variable, 
and at the same time, the regression coefficient vector 
6 ICYL_LS ± as a vector of the regression coefficient 
AICYL__LS± and the constant term BICYL_LSi is calculated 
by the sequential least-squares method. This makes it 
possible to calculate the learned correction value 
KICYL_LS ± as a value in which the variation in intake 
air amount between the cylinders is properly reflected 
therein, even when the engine 3 is in an operating 
condition which can change drastically, such as a 
transient operating condition, causing a sudden change 
in the state of variation in air-fuel ratio between the 
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cylinders. Therefore, even when the engine 3 is in a 
transient operating condition, it is possible to 
properly control the air-fuel ratio while compensating 
for the variation in intake air amount between the 
cylinders . 

Further, the intake air amount variation 
coefficient W ± and regression coefficient vector 6 
ICYL_LSi are calculated by the sequential least -squares 
method, and therefore, differently from the case of 
using the general least -squares method as the 
statistical algorithm, it is possible to calculate the 
intake air amount variation correction coefficient 
KICYLi and the learned correction value KICYL__LSi every 
control cycle, even at the start of the air-fuel ratio 
control. Therefore, by setting the initial values of 
the KICYVi, and KICYV_LSi in advance, it is possible to 
calculate the final fuel injection amount TOUT ± as a 
value always corrected by the product of the intake air 
amount variation correction coefficient KICYLi and the 
learned correction value KICYL_LS ± calculated every 
control cycle, at the start of the air-fuel ratio 
control, whereby the controllability at the start of 
the air-fuel ratio control can be enhanced. This makes 
it possible to enhance the emission reduction rate of 
the catalyst at the start of the air-fuel ratio control. 

It should be noted that in the air-fuel ratio 
control by the air-fuel ratio controller, the intake 
air amount GAIR has characteristics that the amount of 
change in the absolute value thereof can be much larger 
than that of change in the detected air-fuel ratio KACT, 
and in such a case, the amount of change in the 
identified value of the vector t/> (k) of the intake air 
amount variation coefficient identified with the 
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identification algorithm by the equations (60) to (66) 
becomes so large that the control system can be 
unstable. To avoid this, it is only required to 
configure the adaptive observer 61 as shown in FIG. 37. 
That is, it is only required that a filter 61 j (filter 
means) comprised of at least one band-pass filters 1 to 
3 shown in FIG. 39 is provided in the adaptive observer 
61, and the air-fuel ratio control is carried out using 
the filtered value GAIR_F(k) obtained by filtering the 
intake air amount GAIR(k) by the filter 61 j. 

This filter 61j is represented by an equation 
(79) in FIG. 40. In the equation (79), the symbols m* 
and n* represent respective predetermined integers . 
Further, the identification algorithm with which the 
adaptive observer 61 identifies the vector i/> (k) of the 
intake air amount variation coefficient is expressed by 
equations (80) to (86) in FIG. 40. Due to this 
configuration, even when the engine is in an operating 
condition in which the intake air amount GAIR(k) 
changes largely, it is possible to generate the 
filtered, value GAIR_F(k) as a value with reduced range 
of variation in the intake air amount GAIR(k) while 
preserving the necessary information for the 
identification of the vector t/> (k) of the intake air 
amount variation coefficient. Therefore, by 
identifying the vector (k) of the intake air amount 
variation coefficient according to the filtered value 
GAIR_F(k) , it is possible to suppress delay in the 
identification and enhance the accuracy of the same, 
thereby further enhance the stability and response of 
the air-fuel ratio control. 

Although in the second embodiment, the first- 
degree equation is used as the regression equation used 
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in the calculation of the learned correction value 
KICYL_LSi, this is not limitative, but an n-th-degree 
equation (n is an integer equal to or larger than 2) 
may be used. In such a case as well, by calculating 
the regression coefficients and the constant terms of 
the n-th-degree equation by the sequential least- 
squares method, it is possible to obtain the same 
advantageous effects as provided by the second 
embodiment. Further, the learned correction value 
KICYL_LS ± may be calculated by using predetermined 
values set for each of a plurality of operating regions 
in advance, as the regression coefficient and the 
constant term of the regression equation. This can 
reduce the time for computing the learned correction 
value KICYL_LSi, thereby reducing the computational 
load on the ECU 2 . 

Further, although in the second embodiment, the 
I-PD control algorithm is employed as the control 
algorithm for causing the intake air amount variation 
coefficient W ± to converge to the moving average value 
^ave thereof, this is not limitative, but it goes 
without saying that other suitable algorithm may be 
employed. For example, instead of the I-PD control 
algorithm, an IP-D control algorithm (differential- 
preceding PID control algorithm) expressed by equations 
(87) and (88) in FIG. 41 may be employed to calculate 
the intake air amount variation correction coefficient 
KICYLi or a response-specified control algorithm 
(sliding mode control algorithm or back- stepping 
control algorithm) expressed by equations (89) to (91) 
in FIG. 41 may be employed to calculate the intake air 
amount variation correction coefficient KICYLi. Even 
when one of these control algorithms is employed, 
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similarly to the case of using the I-PD control 
algorithm according to the present embodiment , it is 
possible to calculate the intake air amount variation 
correction coefficient KICYL ± such that overshooting is 
not caused in the behavior of the intake air amount 
variation coefficient W ± being caused to converge to 
the moving average value ?ave. As a result, it is 
possible to promptly and appropriately correct the 
variation in intake air amount between the cylinders. 

Further, even when the I-PD control algorithm or 
the response-specified control algorithm is employed in 
the calculation of the intake air amount variation 
correction coefficient KICYLi, as described above, by 
properly setting the feedback gains and the switching 
function-setting parameter to respective values such 
that the speed at which the intake air amount variation 

coefficient ^iCk) converges to the moving average 
value?ave is higher than the speed at which the air- 
fuel ratio variation coefficient Oi(k) converges to 
the moving average value <E>ave, it is possible to 
prevent the air -fuel ratio control by the third air- 
fuel ratio controller 60 and that by the first air-fuel 
ratio controller 30 from interfering with each other. 
Further, each feedback gain thereof may be determined 
based on the optimal regulator theory or the Hoo control 
theory. This makes it possible to more effectively 
suppress overshooting in the behavior of the intake air 
amount variation coefficient W t being caused to 
converge to the moving average value ?ave thereof, 
with the result that the accuracy of correction of 
variation in intake air amount between the cylinders 
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can be further enhanced. 

Further, it goes without saying that when the 
setting time over which the intake air amount variation 
coefficient ? A converges to the moving average value 
ave may be long, the intake air amount variation 
correction coefficient KICYL ± may be calculated with 
the PID control algorithm described above. Further, 
the average value of the intake air amount variation 
coefficient as a target value to which the intake air 
amount variation coefficient ?i is caused to converge 
is not limited to the moving average value ?ave in the 
present embodiment, but it may be a weighted average 
value . 

Further, in the illustrated example of the second 
embodiment, the adaptive observer 61 of the third air- 
fuel ratio controller 60 identifies the vector tp (k) of 
the intake air amount variation coefficient, by the 
variable-gain sequential least-squares method expressed 
by the equations (60) to (65) in FIG. 27, it goes 
without saying that the identification algorithm with 
which the adaptive observer 61 identifies the vector ii> 
(k) of the intake air amount variation coefficient is 
not limited to this. For example, the vector ^(k) of 
the intake air amount variation coefficient may be 
identified by the fixed gain method to which is applied 
the 5 correction method expressed by equations (92) to 
(99) in FIG. 42. 

The symbol i/>base in the equation (92) in FIG. 42 
represents a reference value vector (reference value) 
defined by the equation (93), and four elements of this 
vector, i.e. reference values Vbasel to ?base4 are 
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calculated by searching a table shown in FIG. 4 3 
according to the exhaust gas volume ESV. Further, the 
symbol d^ (k) in the equation (92) represents a 
correction term (correction component) defined by the 
equation (94), and calculated by equations (95) to (99). 

When the vector i/> (k) of the air-fuel ratio 
variation coefficient is identified by the fixed gain 
method to which the 6 correction method is applied, 
the computing time can be reduced compared with the 
case of using the sequential least-squares method, and 
the computational load on the ECU 2 can be reduced. As 
a result , it is possible to reduce the size and cost of 
the ECU 2. Moreover, even when the engine 3 is in an 
operating condition in which the air-fuel ratio changes 
violently, such as a transient operating condition, the 
identified value of the vector ii> (k) can be constrained 
to values close to a value of 1, which makes it 
possible to promptly and properly calculate the vector 
^(k) of the air-fuel ratio variation coefficient 
representative of the variation in intake air amount 
between the cylinders as a value in which the behavior 
of the air-fuel ratio is properly reflected, whereby 
the stability of the air-fuel ratio control can be 
enhanced. 

Although in the preferred embodiments described 
above, the present invention is applied to the control 
system of the engine 3 for automotive vehicles, this is 
not limitative, but it goes without saying that the 
control system according to the present invention can 
be applied to shipborne internal combustion engines and 
other industrial machines . 

It is further understood by those skilled in the 
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art that the foregoing is a preferred embodiment of the 
present invention, and that various changes and 
modifications may be made without departing from the 
spirit and scope thereof. 



